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PREFACE 

The  work  described  in  this  final  report  was  supported  by  the  Advanced  Re¬ 
search  Projects  Agency  of  the  Department  of  Defense  and  was  monitored  by  the 
U.S.  Dureau  of  Mines  under  Contract  No.  H0220064,  ARPA  Order  No.  1579,  Amend¬ 
ment  3  (Program  Code  2F10).  Inclusive  dates  for  this  reporting  period  are  5  May 
1972  through  30  April  1973.  The  Project  Officer  was  Frank  Ruskey.  The  Principal 
Investigator  was  Philip  L.  Jackson.  The  Environmental  Research  Institute  of 
Michigan  number  for  this  report  is  191700 -8 -F. 

The  scientists  who  contributed  to  this  project  throughout  the  contract  period 
are  Ben  Drake,  Robert  Vincent  and  Thomas  Wagner.  Vincent  was  primarily  respon¬ 
sible  for  developing  the  ratio  technique  for  geologic  and  lithologic  problems,  Wagner 
for  soils  analysis,  and  Drake  for  the  geologic  background  and  interpretation.  Jackson 
was  responsible  for  the  overall  coordination  of  the  project.  The  work  reported 
commenced  in  April  1971,  under  Contract  No.H0220041,  ARPA  Order  No.  1579, 
Amendment  2,  Program  Code  IF  10.  This  earlier  work  was  reported  in  Willow  Run 
Laboratories  Report  100 18 -13 -F. 
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GEOLOGIC  RECONNAISSANCE  AND  LITHOLOGIC  IDENTIFICATION 
BY  REMOTE  SENSING 

FINAL  REPORT 

5  May  1972  Through  30  April  1973 


INTRODUCTION 

The  work  during  this  period  was  primarily  concerned  with  a  test  site  for  geologic  remote 
sensing.  After  the  initial  phase,  we  determined  that  the  most  effective  contribution  would  be  a 
test  case.  The  purpose  of  the  test  case  was  to  evaluate  the  usefulness  of  certain  types  of  re¬ 
mote  sensing  and  data  analysis  as  aids  in  determining  the  selection  of  tunnel  sites.  In  effect, 
the  contribution  of  remote  sensing  to  tunnel  site  selection  is  in  our  increasing  knowledge  of 
geology,  lithology,  and  hydrology  (as  described  in  the  previous  report).  The  determination  of 
pertinent  hydrology  by  remote  sensing  can  now  be  considered  operational  (as  illustrated  in 
Section  3.6. 1.1  of  the  previous  report).  Geology  (in  terms  of  tectonics,  depositional  crosional 
patterns  and  ensuing  inferences)  is  the  subject  of  much  literature  in  photogeology  and  side¬ 
looking  radar. 


Lithology,  however,  has  been  the  weakest  area  of  remote  sensing.  Heretofore,  there  has 
been  no  method  of  disciminating  between  rock  types,  except  for  some  research  that  has  been 
done  with  low-flying,  non-imaging  gamma  ray  sensors.  The  power  of  point -by -point  image 
ratioing  of  multispectral  channels  for  both  Si02  and  Fo^  content  of  rocks  was  shown  in  sec¬ 
tion  4.1  of  the  previous  report. 

The  identification  of  rock  types  by  remote  sensing  meant  the  breaking  of  new  ground  aid 
the  filling  in  of  a  gap  in  geologic  remote  sensing.  We,  therefore,  felt  that  our  most  useful  con¬ 
tributions  could  be  made  by  extending  such  identification  and  by  testing  in  an  area  which  in¬ 
cludes  many  types  of  geologic  structures  and  lithologies,  some  of  them  covered  by  vegetation. 

A  test  site  was  chosen  in  the  Black  Hills  of  South  Dakota,  and  airborne  multispectral  scanning 
data  were  obtained  in  May,  1972.  Croundtruthing  at  the  time  of  the  flight  and  after  preliminary 
processing  of  the  data  indicated  a  significant  advance  in  geologic  remote  sensing. 

By  ratioing  pairs  of  multispectral  images,  the  following  geologic  features  can  be  dis¬ 
criminated: 


(1)  exposed  rock  areas 

(2)  separate  rock  groups,  including  Precambrian  quartzite;  Pre-ambrian  phylite,  slate, 
and  schist;  basal  quartzite  and  sandsto.  e:  and  intrusives  such  as  rhyolite,  quartz 
latite,  quartz  monzonite,  and  phonolite 
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(3)  lai  go-scale  structures:  faults,  folds,  intrusive  l)odi(‘s,  and  regional  trends. 

Thu  form  of  the  analysis  output  is  photo-interpretive,  the  cost  is  relatively  low,  the  tech¬ 
nique  has  operational  potential  because  it  is  relatively  simple,  and  no  a  priori  information  is  re¬ 
quired 

To  aid  in  making  the  data  proc  ssme  automatic,  a  gating  logic  has  been  developed  for  the1 
1  utioed  images.  The  basis  of  the  lope  is  to  use  a  known  ratio  range  for  a  rock  type  in  the 
prepai  ation  of  recognition  maps.  All  ratio  values  which  fall  into  a  preset  range  indicating  a 
given  lock  type  are  selected  and  automatically  placed  in  the  location  of  the  rock's  occurrence 
in  'In  scene. 

Statistical  procedures  have  been  investigated  to  indicate  the  probabilities  and  standard  de¬ 
viations  of  rock  type  selections  by  automatic  gating  logic.  The  techniques  are  based  on  labora- 
*'1i  y  measurements  which  enable  statistical  determinations  and  give  a  firm  basis  for  remote 
sensing  of  rock  types. 

The  ratioings  for  SiC^  find  F0O2  ^utent  have  already  been  described  in  t lie  previous  report. 
Much  of  the  present  report  is  concerned  with  establishing  theoretical,  experimental,  and  data 
processing  bases  for  the  ratioing  technique. 
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GENERAL  CONSIDERATIONS  IN  REMOTE  SENSING  FOR  TUNNEL 

SITE  SELECTION 

The  rationale  behind  remote  sensing  for  tunnel  site  selection  is  presented  in  some  detail 
in  the  previous  report.  Also  treated  are  the  aspects  of  geology,  lithology,  and  hydrology  which 
are  required  for  tunned  site  selection  and  for  which  remote  sensing  is  useful. 

During  the  pas  two  years,  remote  sensing  has  advanced  technically  while,  through  wide 
publicity,  the  pot<  ntial  user  has  become  aware  of  its  possibilities.  These  advances  in  techni¬ 
cal  and  information  dissemination  have  come  about  as  a  result  of  the  successful  experiments 
of  ERTS-1. 

The  authors  feel  that  it  is  no  longer  necessary  to  describe  remote  sensing,  point  out  appli¬ 
cations,  and  justify  its  use  for  geologic  and  site* -planning  purposes.  Doth  the  public  and  private 
sector  have  access  to  ERTS  imagery  at  nominal  cost.  And  though  the  fact  is  less  well  known, 
the  public  sector  has  access  to  Strategic  Air  Command  synthetic -aperture  radar  imagery  for 
much  of  the  United  States — through  Goodyear  Aerospace  in  Litchfield,  Arizona.  Also,  applica¬ 
tions  committees  for  synthetic -aperture  radar  are  being  formed  by  NASA  for  EO_  (Earth  Ob¬ 
servation  Satellite)  and,  for  the  Space  Shuttle,  studies  of  multispectral  radar  have  been  made. 

In  the  near  future,  therefore,  photographic,  multispectral,  and  radar  imagery  will  he  avail¬ 
able  to  both  public  and  private  sectors  at  nominal  cost.  Moreover,  since  the  resolution  of 
synthetic -aperture  radar  imagery  is  not  a  function  of  distance,  satellite  imagery  can  he  as  fine 
in  resolution  as  airborne  imagery. 

It  is  likely  that  the  agency  or  private  contractor  responsible  for  a  major  construction  proj¬ 
ect  such  as  a  deep  rock  tunnel  would  use  available  remote  sensing  imagery  for  site  selection. 

In  fact,  it  is  the  opinion  of  the  authors  that  within  a  few  years  specification  forpreconstruction 
geology  studies  will  require  the  full  use  of  remote  sensing  imagery  and  analysis. 

The  cost  of  remote  sensing  is  rapidly  decreasing  for  large  construction  projects.  As  an 

2 

example,  the  Black  Hills  test  case  flown  over  an  80  mi  area  in  South  Dakota  by  the  ERIM  MSS 

2 

aircraft  hased  in  Ann  Arbor,  Michigan,  cost  $100  per  mi  ,  at  10  ft  resolution.  In  this  case,  tran- 

2 

sit  costs  to  South  Dakota  were  shared;  otherwise  the  cost  would  have  been  $150  per  mi  for  a 
single  flight  at  this  distance  from  the  aircraft  home  base.  The  cost  per  square  mile  is  reduced 
when  the  total  area  is  increased. 

2 

The  commercial  rate  for  synthetic -aperture  radar  imagery  is  $6.00-$7.00/mi  ,  based  on 
large-area  coverage.  Specific  flights  over  small  areas  would  Increase  this  cost  many  times. 

Beciuse  remote  sensing  data  have  become  increasingly  available  in  the  last  two  years,  we 
felt  that  we  could  make  the  greatest  contribution  to  geologic  remote  sensing  by  acquiring 
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multispectral  scanning  data  at  a  test  site.  Cost  restrictions  of  the  contract  prevented  the  fund¬ 
ing  of  a  radar  flight  over  the  same  region. 

So,  rather  than  devoting  our  time  to  another  general  analysis  and  evaluation  of  remote  sens¬ 
ing,  this  experiment  was  undertaken  with  concurrence  of  the  project  monitor.  We  feel  fortunate 
that  we  were  able  to  make  significant  advances  in  the  hitherto  weakest  portion  of  geologic  re¬ 
mote  sensing:  lithologic  identification. 
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GEOLOGIC  OUTLINE  OF  THE  NORTHERN  BLACK  HILLS 

The  Black  Bills  of  western  South  Dakota  and  northeastern  Wyoming  are  an  elliptical  dome 
approximately  12 a  miles  long  and  GO  miles  wide  that  trends  approximately  N  35°  W.  The  dome 
actually  consists  of  two  nearly  flat -topped,  steep-sided,  probably  fault -bounded,  uplifted  blocks 
jointed  along  an  approximately  north-south  lineament  passing  near  the  South  Dakota-Wyoming 
border.  The  eastern  block  trends  approximately  north  and  stands  slightly  higher  than  the 
western  block,  which  trends  northwest  from  the  northwestern  corner  of  the  eastern  block.  The 
higher  parts  of  the  Hills  rise  almost  4000  ft  above  the  surrounding  plains. 

The  Black  Hills  basically  consist  of  a  core  of  Precambrian  igneous  and  metamorphic  rocks 
exposed  over  an  area  GO  miles  by  20  miles.  Around  the  core  are  draped  a  few  thousand  feet 
of  Paleozoic  sedimentary  rocks,  which  dip  off  the  dome  in  all  directions.  The  Red  Valley,  with 
Lower  Mesozoic  sedimentary  rocks  underlying  it,  completely  encircles  the  uplift.  A  Cretaceous 
hogback  ridge  forms  an  outer  encircling  rim.  Early  Tertiary  intrusive  rocks  occur  in  a  belt 
across  the  northern  Blark  Hills. 

The  Black  Hills  have  undergone  two  major  periods  of  deformation:  Precambrian  and  early 
Tertiary,  probably  Eocene.  The  Precambrian  deformation  produced  a  series  of  plunging  isocli¬ 
nal  folds  of  various  sizes.  During  or  after  the  closing  phases  of  the  Precambrian  folding,  there 
were  intrusions  of  granite  bodies  and  dioritc  and  gabbro  dikes.  Although  the  granites  occur 
mainly  in  the  southern  Black  Hills,  they  are  exposed  at  four  localities  in  the  northern  Hills.  The 
dioritc  and  gabbro  dikes  were  metamorphosed  to  amphibolites  and  metagabbros  during  the 
Precambrian  deformation:  however,  Bayley  indicates  that  the  amphibolites  are  metamorphosed 
lava  flows,  at  least  in  the  Galena- Roubaix  district  [lj.  Noble  and  Harder  suggest  that  the  pro¬ 
gressive  metamorphism  of  the  rocks  in  the  Lead  region  is  related  to  an  intrusion  of  unexposed 
Precambrian  igneous  rocks  [2 ] . 

The  major  structural  features  of  the  Black  Hills,  as  well  as  the  present  domal  form  of  the 
Hills,  were  formed  in  the  early  Tertiary  during  the  Laramide  Revolution.  Whether  the  Hills 
were  uplifted  by  vertically  acting  forces,  by  laterally  compressive  forres.or  by  some  combina¬ 
tion  of  the  two  types  is  an  open  question;  most  geologists  favor  vertically  acting  forces.  It  has 


1.  Richard  W.  Bayley,  Structure  and  Mineralization  of  Precambrian  Rorks  in  the  Galcna- 
Roubaix  District,  Black  Hills,  South  Dakota,  USGS  Bull.  1312-E,  1970.  15pp. 

2.  J.  A.  Noble  and  J.  O.  Harder,  Stratigraphy  and  Metamorphism  in  a  Part  of  the  Northern 
Black  Hills  and  the  Homestake  Mine,  Lead,  South  Dakota,  Geol.  Soc.  Amer.  Bull.  V.  59,  1948 
pp.  941-976. 
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Ihoii  postulated  that  a  Tertiary  batholith  underlies  the  nortlu  rn  Black  Hills,  if  not  the  entire 
Hills  area,  and  is  responsible  for  the  uplift  [3,  4], 

Several  secondary  structural  features  which  probably  formed  at  the  time  of  the  major  up¬ 
lift  are  superimposed  on  the  general  domal  structure.  These  include  the  following. 

(1)  Sharp  steep  monoclinal  flexures  in  the  sedimentary  rocks,  especially  at  the  western 
edges  of  the  uplifted  blocks  and  along  the  lineament  at  which  the  tv  o  domal  blocks  are  joined  5] . 
The  sedimentary  strata  are  nearly  horizontal  on  top  of  the  dome  (except  where  affected  by  em¬ 
placement  of  the  Tertiary  intrusive  rocks),  but  they  form  monoclines  as  they  dip  off  the  Black 
Hills,  probably  as  a  result  of  draping  over  steep  deep-seated  faults.  In  general,  the  beds  on 

the  western  margins  of  the  blocks  dip  more  steeply  than  those  on  the  (  astern  edges. 

(2)  Small  to  moderate -size  igneous  intrusive  bodies. 

(3)  Plunging  anticlinal  folds,  mainly  at  the  north  ends  of  both  blocks  hut  also  at  the  southern 
end  of  the  eastern  block.  These  plunging  folds  may  be  caused  by  intrusive  bodies  or  by  move¬ 
ments  along  secondary  faults. 

Tertiary  intrusive  igneous  bodies  occur  in  a  r.ont  15  to  20  miles  wide  extending  N  70-75° 

W  across  the  northern  Black  Hills  from  the  Bear  Lodge  Mountains  to  Bear  Butte.  Superim¬ 
posed  on  both  domal  blocks,  these  intrusives  consist  of  relatively  minor  dikes  and  sills,  as  well 
as  larger,  much  more  important  intrusives  (roughly  circular  to  slightly  elongated)  that  range 
from  less  than  a  mile  to  about  10  miles  in  diameter.  It  is  thought  that  the  larger  intrusives  were 
forcibly  injected;  some  consider  them  mainly  laccoliths,  others  describe  them  as  mainly  steep¬ 
sided  stocks  with  concordant  upper  surfaces.  The  larger  intrusives  folded  and  or  faulted  the 
Precambrian-to-Tertiary  rocks  that  they  intruded.  <  renting  a  series  of  domes  that  have  cores 
of  Tertiary  intrusive  rocks.  One  of  the  domes  of  the  northern  Black  Hills  is  the  Lead -Dead  wood 
dome.  It  is  about  12  miles  by  10  miles,  the  longer  dimension  extending  in  a  westerly  direction. 
Near  the  center  of  the  dome  stands  a  very  ragged  stock,  the  Cutting  Stock,  about  2  by  4  miles. 
There  are  also  several  other  smaller  stocks  and  vast  numbers  of  large  and  small  dikes  in  the 
northern  Black  Hills. 


3.  N.  H.  Darton  and  S.  Paige,  Central  Black  Hills  Folio  (219),  I’SGS  Geologic  Atlas,  1925, 

w^pp. 

n,  ,4'„J,',A'  N,Hbk’’  J'  °'  Ilardcr>  and  A-  L-  Slaughter,  Structure  of  a  Part  of  the  Northern 
Black  Hills  and  the  Homestake  Mine,  Lead,  South  Dakota,  Geol.  Soc.  Amer,  Bull.  V.  GO,  1949, 
pp.  321  -352.  1  * 

5.  L.  II.  Shapiro  and  ,J.  P.  Cries,  Ore  Deposits  in  Rocks  of  Paleozoic  and  Tertiary  Ace  of 
kilfs’  South  Dak°ta,  L’SGS  Heavy  Metals  Program,  Contract  No.  14-08-0001  - 

10650,  19/0,  235pp. 
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1!.  W .  Bayley  of  the  L'.S.  Goo  logical  Survoy  lias  just  finished  a  seven-year  mapping  project 
(1965-1972)  in  which  he  revised  the  Prerambrian  geology  of  the  northern  Black  Hills.  The 
results  of  his  mapping  are  available  as  an  open  file  map  and  two  publications  [1,  6].  Except  for 
the  Lead  region,  the  Prccambrian  geology  of  the  northern  Black  Hills  generally  was  not  well 
known  before  Bayley's  mapping.  Bayley  has  made  the  only  overall  correlation  of  Prccambrian 
rocks  in  different  parts  of  the  northern  Hills,  although  some  correlations  between  differenl 
districts  had  been  done  previously. 

Bayley's  recent  mapping  in  the  northern  Black  Hills  has  shown  that  the  Precanibrian  rock 
types  present  include  granite,  graywarke. slate,  phyllile,  schists,  chert,  quartzite,  conglomerate, 
arkose,  limestone,  banded  iron  formation  (taconite),  metagabbro  and  or  amphibolite,  and  meta¬ 
basalt.  Three  well-defined  zones  of  metamorphisni  in  the  Precanibrian  rocks  of  the  Load  region, 
which  have  been  described  by  Noble  and  Harder  f 2 1 ,  have  an  increasing  degree  of  metaniorphism 
from  southwest  to  northeast. 

The  sedimentary  rocks  in  the  northern  Black  Hills  mainly  range  in  age  from  Late  Cambrian 
to  Early  Cretaceous,  but  some  are  as  young  as  Tertiary.  The  Phanerozoic  formations  consist  of 
sandstone  and  quartzite,  shale,  limestone,  dolomite,  siltstone,  days,  conglomerate,  and  some 
thin  coal  seams. 

At  many  places  around  the  eastern  side  of  the  northern  Black  Hills  (and  tit  scattered  local¬ 
ities  in  the  Hills)  are  deposits  of  the  Middle  Oligocene  White  ltiver  Group  laid  down  by  streams 
after  Ihe  uplift  of  Ihe  Black  Hills.  These  deposits  are  i  rosional  remnants  of  a  former  wide¬ 
spread  mantle  of  sand  and  coarse  gravel  underlain  in  places  by  fuller's  earlh  and  clay.  Also 
around  the  northern  Black  llills  are  different  levels  of  Quaternary  sand,  gravel,  and  loam  terrace 
deposits.  Quaternary  alluvium  occurs  in  the  major,  wider  stream  valleys  bolh  in  and  around  Ihese 
northern  Hills. 

Tertiary  acidic  volcanics,  probably  Eocene  in  age,  are  found  in  at  least  four  general  areas 
in  Ihe  northern  Hills. 

No  comprehensive  study  has  been  done  of  the  petrology  and  petrography  of  the  Tertiary 
intrusives  in  the  northern  Black  Hills.  Identified  intrusive  rock  types  include  rhyolite,  quartz 
latite,  monzonite,  quartz  monzonite,  phonolitc,  and  grorudite — porphyries  of  these  intrusive 


6.  Richard  W.  Bayley,  Preliminary  Geologic  Map  of  the  Nemo  District,  Black  Hills, 
South  Dakota,  Misc.  Geologic  Investigations,  Map  1-712,  1972. 
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typos  are  common.  There  seem  to  he  gradations  from  one  type  to  another.  Most  of  the 
larger  intrusives  are  multiple  intrusions  of  several  different  types  of  rock;  and  it  seems  safe 
to  say  that  most  of  the  multiple  intrusion  bodies  have  not  been  mapped  correctly. 

The  structural  geology  of  the  northern  Black  Hills  is  known  in  detail  only  for  scattered 
localities  where  there  has  been  concentrated  mapping,  such  as  Lead,  Rorhford,  the  Nemo  dis¬ 
trict,  and  the  Galena -Roubaix  district.  It  is  known  that  there  is  a  series  of  large,  isoclinal 
anticlines  and  synclines  with  smaller  isoclinal  folds  on  their  limbs  in  the  Precambrian  rocks 
of  the  northern  Black  Hills,  and  that  all  these  folds  were  formed  during  the  Precambrian  de¬ 
formation.  These  large  folds  generally  trend  north  to  northwest  and  their  limbs  dip  steeply 
northeast:  the  folds  plunge  at  angles  that  are  gentle  to  almost  vertical.  Sinuous  refolded  iso¬ 
clinal  folds  (cross  folds)  occur  in  several  regions  of  the  northern  Hills.  The  age  of  the  cross 
folding  is  not  known;  it  might  be  early  Tertiary  instead  of  Precambrian.  Slaughter  believes  that 
tlu  structural  effects  of  the  injection  of  Tertiary  intrusive  rocks  were  very  important  [7|. 

Bayley's  mapping  has  shown  the  presence  of  several  large  generally  northwest -trending, 
right -lateral,  strike-slip  faults.  Major  normal  faults  also  have  been  mapped,  as  well  as  rel¬ 
atively  few  reverse  faults;  these  reverse  faults  generally  merge  into  other  types  of  faults. 


7.  A.  L.  Slaughter,  The  Homcstake  Mine;  Ore  Deposits  in  the  United  States  1933-1967 
AIME,  1968,  pp.  1436-1459. 
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RATIO  GATING  LOGIC  RECOGNITION  MAPS  PRODUCED  I  ROM  ANALOG 

RATIO  IMAGES 

Ratio  images  employed  for  the  first  time  two  years  ago  under  a  NASA  task  have  been  found 
usetul  lor  geolog  .d  remote  sensing  [8 -12} .  Ratio  images  are  uselul  because  they  (1)  enhance 
the  contrast  between  targets  which  have  intersecting  spectral  reflectance  curves:  (2)  suppress 
or  reduce  illumination  variations  across  the  scene;  and  (3)  displa.  spectral  data  from  all  points 
in  the  scene  before  automatic  classifications  have  been  made  that  result  in  some  areas  being 
classified  unrecognized.  In  addition,  they  are  relatively  simple  for  the  photointerpreter  to 
use  —because  of  their  similar  appearance  to  aerial  photos  and  the  relative  ease  by  which  they 
can  be  qualitatively  related  to  laboratory  spectra.  Ratio  images  also  maintain  the  advantages 
of  niultispectral  scanner  in  that  they  have  greater  spectral  range,  are  easier  to  process 
quantitatively,  and  require  less  mosaicking  than  aerial  photos. 

There  are  still  needs,  however,  for  recognition  maps  in  geological  remote  sensing.  Orig¬ 
inally,  the  ratio  gating  logic  described  in  this  section  was  created  simply  to  enable  the  geologi¬ 
cal  user  to  produce  recognition  maps  of  some  targets  without  having  to  forfeit  the  processing 
time  already  invested  in  producing  the  ratio  images.  Traditional  analog  niultispectral  recogni¬ 
tion  methods  utilizing  likelihood -ratio,  statistical,  decision-theory  require  several  things  be¬ 
fore  recognition  maps  are  produced;  a  description  of  training  sets  (including  location  and 
spectral  characteristics)  by  the  user  to  a  SPARC  operator;  the  preparation  of  magnetic  tape 
loops;  and  training  of  the  SPARC  computer  on  the  training  sets.  In  short,  the  user  of  ratio 
images  would  have  to  star  from  scratch  to  produce  analog  recognition  maps  in  this  manner. 

In  the  course  of  this  investigation,  however,  several  advantages  of  the  ratio  gating  logic  over 
the  traditional  niultispectral  recognition  method  have  become  more  obvious,  so  this  rather 


8.  It.  K.  Vincent,  R.  Horvath,  F.  Thomson,  and  E.  Work,  Remote  Sensing  Data-Analvsis 
Projects  Associated  with  the  NASA  Earth  Resources  Spectral  Information  System,  The  Uni¬ 
versity  of  Michigan,  Report  No.  3 165 -77 -T,  NASA  Contract  NAS9-9784,  1972. 

9.  R.  K.  Vincent  and  F.  Thomson,  Spectral  Compositional  Imaging  of  Silicate  Rocks,  J.  of 
Geophys.  Res.,  Vol.  77,  1972,  pp.  2465-2471. 

10.  R.  K.  Vincent,  Rock-Type  Discrimination  from  Ratio  Images  of  the  Pisgali  Crater, 
California  Test  Site,  The  University  of  Michigan,  Report  No.  3165 -77 -T,  NASA  Contract 
NAS9-9784,  1972. 

11.  It.  K.  Vincent,  F.  Thomson,  and  K.  Watson, Recognition  of  Exposed  Quartz  Sand  and 
Sandstone  by  Two-Channel  Infrared  Imagery,  J.  of  Geophys.  Res.,  Vol.  77,  1972,  pp.  2473-2477. 

12.  R.  K.  Vincent  and  F.  Thomson,  Rock  Type  Discrimination  from  Ratioed  Infrared  Scan¬ 
ner  Images  of  Pisgali  Crater,  California,  Science,  Vol.  175,  1972,  pp.  986-988. 
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simpk  technique  has  taken  on  greater  significance.  Although  it  theoretically  rannot  be  as  power¬ 
ful  an  automatic  derision-makinR  method  as  the  more  complex  statistical  ones,  its  advantages 
out wei rIi  the  costs  in  decision  making  for  those  remote  sensing  problems  characterized  by  rel¬ 
atively  hiRh  spectral  contrasts  either  ainonK  several  targets  of  interest  or  between  tarots  and 
background. 

In  Section  4.1,  the  recognition  procedure  involving  ratio  gating  logic  (hereinafter  called 
RAGAL)  will  be  described,  followed  in  Section  4.2  by  the  experimental  results  from  its  first 
application  to  scanner  data.  In  Section  4.3,  traditional  multispectral  recognition  results  will 
be  Riven  and  compared  with  RAGAL  results.  Finally,  in  the  discussion  and  conclusions  sec¬ 
tion  (4.4),  the  advantages  and  disadvantages  of  the  RAGAL  procedure  will  be  listed,  alo.iR  with 
suggestions  for  future  improvement. 


4.1  DESCRIPTION  OF  TIIE  RAGAL  PROCEDURE 

Analog  ratio  images  are  produced  by  inputting  signals  from  two  channels  of  a  multispec¬ 
tral  scanner  into  an  analog  computer,  which  first  modifies  the  signal  voltage  in  each  channel 
to  be  directly  proportional  to  detected  radiance  and  then  divides  the  voltage  of  one  channel  by 
that  of  the  other.  The  resulting  ratio  voltage,  which  is  directly  proportional  to  the  ratio  of  ra¬ 
diances  detected  in  the  two  channels,  is  then  fed  to  a  glow  tube  that  scans  over  a  moving  strip 
of  photographic  film.  Film  density  on  the  resulting  photographic  print  is  a  direct  function 
oi  the  ratio  of  radiances  in  the  two  channels  of  interest.  A  grayscale  is  printed  with  each 
ratio  image  to  enable  the  user  to  relate  film  density  to  voltage,  which  is  proportional  to  the 
ratio.  Whenever  ratio  images  are  produced,  the  electronic  settings  required  on  the  SPARC 
ratio  processor  are  recorded,  so  that  their  reproduction  at  a  later  time  is  possible. 

The  user  can  treat  ratio  images  as  an  end  product  in  themselves.  The  RAGAL  method  is 
designed,  however,  to  permit  the  user  to  make  recognition  maps  based  on  the  ratio  images,  if 
he  so  desires.  If  the  user  defines  "m"  training  sets  and  has  "n"  ratio  images,  he  must  deter¬ 
mine  'n  voltage  ranges  for  each  of  the  "m”  training  sets.  A  training  set  ran  consist  of  any 
area  m  the  ratio  image  or.  if  RAGAL  can  be  developed  to  its  final  intended  form,  of  a  group 

of  laboratory  spectra.  At  the  present  time,  only  the  former  type  of  training  set  has  been  im¬ 
plemented. 

Once  a  target  of  interest  has  been  defined  in  the  ratio  imagery,  the  user  determines  the 
target  voltage  range  for  each  of  the  "n"  ratio  images  by  comparing  the  film  density  over  the 
target  area  with  the  corresponding  grayscale  for  each  ratio.  This  can  be  accomplished  either 
by  visual  examination  or,  more  accurately,  by  a  photodensitometer.  The  resulting  set  of  volt¬ 
age  ranges  for  each  target  is  communicated  to  the  SPARC  operator,  who  first  reproduces  the 
original  electronic  settings  used  to  construct  the  ratio  images  and  then  creates  a  series  of 
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A\0  gates  on  the  SPARC.  Thorn  is  an  AND  gate  for  each  of  the  "n"  ratios,  with  the  upper  and 

lower  voltage  (proportional  to  radiance  ratio)  values  on  each  gate  defined  by  the  voltage  ranges 

ol  a  givm  target.  For  any  point  in  the  scene  to  he  recognized  as  that  target,  the  ratios  of  that 

point  must  fall  within  the  limits  of  every  one  ol  the  AND  gate's.  All  points  meeting  these  criteria 

an  printed  as  white,  and  all  non -target  points  are  printed  black.  Tlu  AND  gates  are  reset  for 

each  different  target.  Since  statistical  likelihood  ratio  decision  rules  are  not  used,  there  can  he 

legions  of  taige't  overlap  an  area  in  the1  scene  can  be*  recognized  as  more*  than  one'  target. 

However,  with  RAGAL  there  is  no  limit  to  the  number  of  targets  that  can  be  attempted  because 

each  target  is  independent.  On  the  other  hand,  SPARC  can  implement  the  statistical  likelihood 

4  8 

ratio  decision  rule  only  for  a  maximum  of  —  interdependent  targets,  where  "n"  is  the  number 
of  input  channels. 

It  was  not  possible  to  implement  RAGAL  exactly  as  desired  on  the  SPARC  computer  be¬ 
cause.  in  its  present  configuration,  SPARC  is  not  capable  of  simple  gating.  Therefore,  a  Euclid¬ 
ian  distance  criterion  was  used  which  approximates  the  gating  logic.  For  a  given  mean  voltage 
and  given  variance  in  a  single  input  channel,  the  Euclidian  distance  rule  weights  the  voltage 
from  a  given  point  in  the  scene  according  to  its  deviation  from  the  mean  voltage.  Therefore,  an 
approximation  to  the  ratio  gate  voltages  was  mane  by  the  mean  voltage  (the  average  of  the  upper 
and  lower  ratio  voltage  limits)  and  a  variance  that  was  1.5  times  the  difference  between  the 
mean  and  lower  (or  upper)  limit  voltages.  Figure  1  is  a  schematic  example  of  this  approxima¬ 
tion. 


4.2  EXPERIMENTAL  RESULTS 

The  preliminary  testing  of  RAGAL  has  been  with  a  data  set  from  a  region  northwest  of 
Rapid  City,  in  the  northern  Black  Hills  of  South  Dakota.  Data  from  this  region  were  collected  by 
the  ERIM  (  -47  aircraft  for  the  U.S.  Bureau  of  Mines  on  May  21,  1972  beginning  at  approximately 
9  A.M  Mountain  Daylight  Time.  Altitude  was  5000  ft  above  terrain.  The  sky  was  tree  of  clouds, 
but  a  low  haze  was  prose  .t  throughout  the  flightline,  which  was  about  24  mi  in  length  and  2  mi  (see 
Fig.  2)  in  width.  Only  visible  and  reflective-IR  scanner  channels  were  ratioed  because  of  a 
malfunction  in  the  two-element  thermal  detector  used  for  '"0Si02  and  "MSiOg-AI^)  discrimina¬ 
tion  (see  Section  5  of  this  report). 

Figures  3  and  4  show  six  ratio  images  of  flightline  section  2A  over  Lead,  South  Dakota. 
Figures  5  and  6  show  the  same  ratio  images  for  section  2B,  the  second  half  of  that  flightline. 
North  is  approximately  to  the  upper  left  corner  in  each  figure.  The  six  spectral  channels  used 
for  these  ratios  are  given  in  Table  1,  where  Rij  represents  the  radiance  in  the  i-th  channel 
divided  by  the  radiance  in  the  j-th  channel.  Considerable  information  is  available  from  the 
ratio  images  alone.  However,  in  this  report  the  topic  of  interest  is  computer  target  recognition, 
which  in  the  future  may  become  a  truly  automatic  function. 
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Vrim 


Variance  =  1.5(Vu  -  Vm) 


+  1 .5 (V  -  V  ) 
u  m 

=  Desired  Upper  Voltage  Limit 
=  Mean  Voltage 

=  Desired  Lower  Voltage  Limit 

-  1 .5 ( V  -  V  ) 
u  m 


FIGURE  1.  EUCLIDIAN  DISTANCE  APPROXIMATION  TO  RATIO  GATING  LOGIC 
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FIGURE  2.  LOCATION  OF  MSS  FLIGHT  LINE  FLOWN  IN  NORTHERN  BLACK  HILLS 
SOUTH  DAKOTA,  ON  MAY  21,  1972.  Sections  2A  and2B  of  the  Hlghtline  are  both  shown. 
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FIGURE  4.  ANALOG  IMAGES  OF  RATIOS  R45.  R^.  AND  R  FOR  FLIGHT  SECTION  2A,  NEAR  LEAD.  SOUTH  DAKOTA 
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FIGURE  5.  ANALOG  IMAGES  OF  RATIOS  R„„,  R  AND  R  FOR  FLIGHT  SECTION  2B,  NEAR  LEAD,  SOUTH  DAKOTA 


25 


FIGURE  R  ANALOG  IMAGES  OF  RATIOS  R  R  ,,  AND  R. .  F^R  FLIGHT  SECTION  2B,  NEAR  LEAD.  SOUTH  DAKOTA 
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TABLE  1.  LIS1  OF  SPECTRAL  CHANNELS  AND  RATIOS  USED  WITH  THE  BLACK  HILLS  DATA 


Single  Channel  No. 

Wavelength  Region  (um) 

Ratios  Including  that  Channel 

1 

.50-. 54 

R41 

2 

.62-. 70 

R32 ’  R42’  R62 

3 

.67-. 94 

R32 

4 

1.0-1. 4 

R41*  R42 ’  R45 

5 

1. 5-1.8 

R45’  R65 

6 

2. 0-2. 6 

R62 ’  R65 
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Figure-  i  is  a  flow  chart  describing  the  attempted  final  categorization  of  targets  in  the  Black 
Bills  data.  The  extremum  ratio  ranges  for  each  of  these  targets  are  given  in  Table  2  On  the 
nasis  of  the  listed  number  ol  test  areas  given  in  the  last  column,  these  are  the  target  c  lassifi - 
cations  which  appear  to  be  uniquely  diseriminable  from  one  another.  Other  target  classifications 
wc  re  examined,  but  their  ratio  ranges  were  found  to  overlap  the  ranges  of  other  targets.  For 
instance,  an  effort  was  made  to  subcategorize  the  igneous  intrusives  into  rhyolite,  quartz  lutite, 
and  phonolite;  but  they  overlapped  each  other  enough  to  force  them  into  a  single,  more  general 
category. 

BAGAL  recognition  maps  of  the  10  targets  listed  in  Table  2  are  shown  in  Figs.  8  through 
10  for  flight  line  section  2A;  for  section  2B,  they  are  shown  in  Figs.  11  through  13  The  top¬ 
most  image  strips  in  Figs.  8  and  11  are  Ii32  ratio  images  of  sections  2A  and  2B,  respectively, 
upon  which  the  test  sites  enumerated  in  Table  2  have  been  indicated.  For  instance,  there  arc 
five  vigorous -range -grass  targets  labeled  ”1."  five  quartzite  targets  labeled  "7,”  etc.  on  see- 
.ions  2A  and  2R  combined.  Target  5,  all  non -vegetation,  was  a  simple  one -ratio  threshold 
Process  for  which  no  specific  training  set  was  used.  That  R32  threshold.  2.2  volts,  was  selected 
on  the  basis  of  data  in  the  first  column  of  Table  2,  in  which  a  natural  separation  occurs  between 
vegetative  and  non-vegetative  targets  at  2.2  volts. 

Three  types  of  experimental  artifacts  must  be-  pointed  out.  First,  the  dark  vertical  regions 
in  all  of  the  section-211  recognition  maps  and  ratio  images  including  the  ratios  R^,  h  ( 

Rgg,  and  R^  were  caused  by  a  temporary  failure  of  the  three-element  detector  responsible  for 
spectral  channels  4,  5,  and  G.  Second,  in  ratio  images  and  vigorous -range -grass  (target  1) 
t  ecognition  maps  of  both  flightline  sections,  some  areas  of  the  vigorous  range  grass  are  brighter 
than  the  brightest  grayscale  level  of  ratio  R^.  This  is  typified  in  R4J  ratio  images  by  black 
regions  and  in  the  target  1  recognition  map  as  unidentified  (black)  vigorous -range-grass  test 
sites.  Only  the  most  vigorous  grasses  were  unrecognized  in  this  manner,  however.  This  problem 
is  trivial  and  can  be  corrected  in  future  applications  of  RAGAL.  Third,  there  is  some  noise  in 
the  target -1  recognition  map  of  section  2A  (white  vertical  streaks)  caused  by  non-recurring 
tape  recorder  noise  while  the  recognition  map  was  being  generated. 

In  addition  to  the  above  errors,  the  reader  is  again  reminded  that  the  human  eye  was  used  to 
determine  voltages  from  the  ratio  image  gray -scales,  and  that  the  gating  logic  was  approximated 
by  a  Euclidian  distance  method.  The  former,  which  was  probably  accurate  only  to  approximately 
+  0.2  volt  (with  a  bias  toward  wider  ratio  ranges  so  the  target  cannot  be  missed),  can  be  greatly 
improved  with  a  photodensitometer.  The  latter  error,  which  is  impossible  to  estimate,  can  In- 
corrected  by  building  a  cheap  AND-gate  processor  to  replace  the  SPARC  logic  for  RAGAL 
application. 

A  brief  evaluation  of  each  of  the  10  target  recognition  maps  is  given  in  Table  3.  The  worst 
recognition  occurred  for  targets  1,  4,  and  10.  The  target-1,  vigorous -range -grass  errors  arc- 
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TABLE  2.  RATIO  RANGES  OF  BLACK  HILLS  TARGETS 


Tar gut 
No. 

Target  Name 

Ratio 

■t - - 

Voltage 

Limits  (Volts) 

No.  of 
Tes  t 

R32 

RA2 

R62 

RA5 

R65 

rai 

Areas 

1 

Vigorous  Range 
Grass 

6. 6-7.0 

A. 9-5. 3 

3.7-A.3 

A.1-A.7 

2. 9-3. 5 

5. 3-5. 5 

5 

2 

Dormant  and 

2. 2-3. A 

2. 0-3.0 

2. A-2.8 

2.3-3. 1 

3. 5-A. 3 

3  5-A  3 

7 

Low  Vegetation 

3 

Conifer  Trees 

A. 3-5.9 

2. 7-3. 7 

2.6-3. 6 

3.5-A.5 

2. 7-3.5 

3.6-A.A 

3 

A 

Deciduous  Trees 

3.0-A.2 

2. 6-3.6 

2. 9-3. 9 

2. 3-3.1 

3.A-A.2 

3.7-A.5 

1 

5 

All  Non- 
Vegetation 

0.6-2. 2 

- 

- 

- 

- 

- 

- 

6 

Water 

0.6-2. 2 

0.2-0. 8 

0.5-2. 5 

1. A-2.8 

0.8-3. A 

0.9-2. 3 

A 

7 

Quartzite 

0.6-2. 2 

0.A-2.0 

2. A-2.8 

1. 3-2.1 

3.7-A.9 

0. 7-1.9 

5 

8 

Slate,  Phyllite, 
Schist 

0.6-1.  A 

0.5-1. 3 

0.A-1.A 

1.5-2. 7 

3. 7-5. 5 

0. 8-3.0 

5 

9 

Igneous  Intrusives 

0.7-1. 3 

1.0-1. 2 

2.0-3. 0 

1. 3-1.9 

A. 7-5. 5 

2.2-2 .6 

7 

(Rhyolite,  Quartz 
Latite,  &  Phonolitc 

) 

10 

L 

Deadwood  Sandstone 

0.6-2. 2 

1. 7-1.8 

2. 2-2. 7 

2.0-2. 5 

A. 3-5. 5 

3. 0-3. 2 

2 
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FIGURE  7.  TARGET  CATEGORIES  FOR  BLACK  HILLS.  SOUTH 
DAKOTA  DATA  SET 
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FIGURE  9.  RATIO  IMAGE  (R32)  AND  RATIO  GATING  LOGIC  RECOGNITION  MAPS  OF  ALL  NON -VEGETATION 
AND  WATER  IN  FLIGHT  SECTION  2A  OF  BLACK  HILLS  DATA 
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FIGURE  11.  RATIO  GATING  LOGIC  RECOGNITION  MAPS  OF  VEGETATIVE  TARGETS  IN  FLIGHT  SECTION  2t  OF 


k 
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FIGURE  12.  RATIO  IMAGE  (R32)  AND  RATIO  GATING  LOGIC  RECOGNITION  MAPS  OF  ALL  NON -VEGETATION 
AND  WATER  IN  FLIGHT  SECTION  2B  OF  BLACK  HILLS  DATA 
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underslood,  however,  and  ran  easily  he  corroded  ,as  explained  earlier).  Errors  target  1 
ead.ood  sandstone)  are  understandable  In  that  the  ratio  ranees  in  11  and  R  for  II  ,  , 

zz  :r  m  ~  °-2  — »■  - 

.2  1,  ol  the  meat,  vollaees)  this  laves, ieation.  Errors  In  taryol  4  are  the  most  dimcnl, 

'  S°  °'  ",C  S"'a"  "mo“  °<  B™nd  'm>"  available  tor  deciduous  trees  and  the 

realistic  overlap  between  deciduous  trees  and  broad-leaf  low  ,eee!a,io„  (part  of  tareet  2,  R 

oenttton  of  both  targets  4  and  10  can  be  considerably  improved  by  the  use  of  a  photodensilomeTer 
and  a  more  exact  tntpl, erne, ttafion  o,  fhe  ra„o  gating  loyle,  without  the  Euclidian  distance  Z!Z- 

,j:rr  7777?? m  **  -  *  «»*** . 

•  k  .  ^32'  In  fact-  Tab,e  2  indicates  that  a  simple  thresholding  0f  ratio  R  may  be 

su  icient  to  discriminate  among  targets  1,  3,  5,  and  a  combined  target  of  2  and  ^Tliis  will 

*  PUr  ‘her  111  ‘hC  n0Xt  RAGAL  investigation.  Other  recognition  maps  winch  appear  to 

b0  good  are  those  for  targets  2,  7,  and  8.  Recognition  maps  for  targets  3  and  6  were  doJn- 

co  7  ,  fa'Se  alarmS’  Wh‘Ch  ^  b°th  CaS6S  WerP  ',ri,11arily  targets  that  had  been  re¬ 

cognized  in  other  recognition  maps. 

4  3  COMPARISON  OF  TRADITIONAL  MULT  (SPECTRAL  RECOGNITION  WITH  RATIO  GATING 

In  order  compare  RAGAL  results  with  tradifiona,  analog  tnultispoc.ra,  recognition  tech- 
ntques,  recognition  maps  of  six  targets  were  produced  for  fliyhtiine  section  2B  w-tlMhe  , 
ttnuca,  iihelihood  ratio  decision  rule.  Maymotic  ,n,»  he,,  ioo.L  fcontaC 7  ,  ^ 

»ere  prepared  for  each  fhe  six  faryets  in  iieu  of  usiny  the  ra„o  in, aye.  The  SPARC  „ 
or  then  recorded  the  input  voitayes  for  each  input  channel  over  each  fes,  site.  Two  sets 
Wttf.on  maps  were  made  in  Ibis  manner,  one  with  single-channel  radiances  as  input,  and 

Z  ”,  'T  ",  ra"°  'nPU'5  <ra"°  '’rCpr"C“",B)'  F,E“™S  14  »*  >5  »bow  the  results  of  usiny 
s  ny  e  c  tannels  of  Table  1  as  inputs.  Fiyuros  16  and  17  show  the  results  of  usiny  ratio 

r::r tor  ,he  ragal  pr”codur"  ^ 4  —  • 

in  „i  „  C°"‘'"0n  S'n"  01  "'B  ’calntne  sites  Rape  loop,,  are 

yhtl.no  section  2A,  It  is  possible  for  a  target  to  be  absent  completely— as  in  the  slate 
phyllite,  schist  case  (target  Thn  t  ^ 

E  *  71  'raining  sets  used  were  common  to  all  three  (vires  of 

the  RAGAL  nlpm  ^  <*"*«  *>•  -re  different  for 

zrrr r  - - 

TTT  ""  C,,eC'  ~  environmental  factors,  such  as  Jhadow '“y  "ZTna 
-on,  and  scan  anyie  variations.  Tile  statistical  dec, si.  theory  with  ratio  prepr.eZiny TI 


FIGURE  14.  RECOGNITION  MAPS  PRODUCED  BY  STATISTICAL  DECISION  THEORY  WITH  SINGLE-CHANNEL  INPUTS 
FOR  VEGETATIVE  TARGETS  IN  FLIGHT  SECTION  2B  OF  BLACK  HILLS  DATA  SET 
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FIGURE  17.  RECOGNITION  MAPS  PRODUCED  BY  STATISTICAL  DECISION  THEORY  WITH  RATIO  INPUTS  FOR 
GEOLOGIC  TARGETS  IN  FLIGHT  SECTION  2B  OF  BLACK  HILLS  DATA  SET 


TABLE  4.  QUALITATIVE  COMPARISON’  OF  RECOGNITION  RESULTS  OF  STATISTICAL  DECISION  THEORY 

WITH  RATIO  GATING  LOGIC 
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best  for  target  0,  but  RAGAL  was  best  for  targets  2  and  H.  Targets  1,  3,  and  7  were  different 
for  these  two  types  of  recognition  logic,  but  neither  type  could  be  called  superior  for  those 
targets.  Therefore,  the  statistical  decision  theory  with  ratio  preprocessing  and  RAGAL  were 
approximately  equal  in  this  experiment.  The  former  probably  could  have  been  further  improved 
by  using  the  same  number  of  test  sites  for  each  target  as  were  used  with  the  latter  method 
However,  that  w  mid  have  required  more  tape  loops  and  SPARC  time  than  could  be  afforded  for 
this  investigation. 

This  investigation  was  too  economically  limited  to  permit  a  more  exhaustive,  quantitative 
comparison  of  these  recognition  methods.  A  better  comparison  can  be  made  in  the  future  when 
the  RAGAL  procedure  is  past  the  prototype  stage. 

4.4  DISCUSSION  AND  CONCLUSIONS 

A  new  method  for  producing  automatic  recognition  maps  has  been  presented  and  applied 
to  scanner  data  for  the  first  time.  The  ratio  gating  logic  (RAGAL)  method  involves  the  fol¬ 
lowing:  (1)  production  of  selected  analog  ratio  images,  (2)  the  user's  determination  of  a 
series  of  ratio  ranges  for  each  target  from  the  analog  ratio  images  and  their  gray-scales, 

(3)  setting  these  ranges  for  each  ratio  on  the  SPARC  computer  to  form  a  series  of  ratio  AND 
gates,  and  (4)  the  production  of  an  automatic  recognition  map  for  each  target  via  AND-gate 
logic.  This  procedure  is  primarily  useful  for  targets  that  have  distinctive  reflective  minima 
(absorption  bands)  related  to  chemical  composition.  Since  one  of  the  most  useful  parameters 
in  geologic  mapping  is  chemical  composition,  this  method  is  particularly  useful  for  geologic 
(and  possibly  soil)  mapping. 

The  procedure  was  not  implemented  exactly  as  desired  in  this  first  application,  but  the 
two  main  sources  of  error  can  be  diminished  in  future  applications  of  RAGAL  by  using  a 
photodensitometer  instead  of  visual  inspection  to  compare  ratio  images  with  grav-scales  and 
by  substituting  a  simple  AND-gate  logic  black  box  for  the  SPARC  computer  to  produce  auto¬ 
matic  recognition  maps  with  extremum  ratio  ranges.  An  alternative  to  the  last  step  would 
he  the  substitution  of  ratio  means  and  standard  deviations  for  extremum  ratio  ranges,  such 
that  the  Euclidian  distance  decision  logic  of  the  SPARC  would  be  more  suitable  than  it  is  now. 

The  RAGAL  procedure  was  qualitatively  compared  with  statistical  decision  theory  (for 
two  types  of  inputs)  by  producing  recognition  maps  of  several  geologic  and  vegetative  targets. 

In  general,  RAGAL  performed  about  as  well  as  statistical  decision  theory  with  ratio  prepro¬ 
cessing,  and  it  was  much  better  than  statistical  decision  theory  with  no  preprocessing  (single- 
channel  radiance  inputs).  The  data  set  used  for  this  comparison  is  from  an  area  in  South 
Dakota  with  sharp  topographic  relief,  which  may  in  part  be  why  ratio  preprocessing  (which 
suppresses  illumination  variations  induced  by  changes  in  slope)  improved  the  automatic  de¬ 
cisions.  In  this  test  the  same  scanner  channels  were  used  for  all  three  types  of  recognition 
maps.  45 
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Although  the  RAGAL  procedure  is  not  as  powerful  a  decision  maker  as  the  statistical 
decision  theory,  it  has  several  advantages  for  the  types  of  targets  sought  in  the  Black  Hills 
data  set’ 

(1)  RAGAL  results  are  relatively  easy  to  relate  to  laboratory  spectra  because  of  the 
ratio  suppression  of  environmental  effects  (slopes,  shadows,  atmospheric  effects, 
etc.)  and  because  of  the  simple  logic  involved. 

(2)  Because  of  this  relatability  to  laboratory  spectra  and  because  the  ratio  images  allow 
the  user  to  determine  the  ratio  values  of  a  target  prior  to  making  an  automatic  rec¬ 
ognition  map,  the  RAGAL  procedure  is  more  predictive  in  the  sense  that  the  user 
has  a  good  idea  about  what  targets  can  and  cannot  be  discriminated  before  the  auto¬ 
matic  recognition  maps  are  produced. 

(3)  Feature  selection  is  relatively  simple  because  the  user  can  determine  what  ratio 
pairs  are  best  for  which  target  from  lab  data,  instead  of  from  scanner  data. 

(4)  Ratio  preprocessing,  an  important  part  of  RAGAL,  should  yield  more  repeatable  re¬ 
sults  from  one  data  set  to  another  than  any  decision  process  based  on  single-channel 
radiances.  This  has  yet  to  lx?  proven,  however 

(5)  Ratio  images,  also  an  important  constituent  of  RAGAL,  are  useful  research  tools  for 
seeking  new  physical  phenomena  to  monitor  by  remote  means.  Patterns  which  the 
observer  can  see  in  ratio  images  are  sometimes  unrecognizable  in  single-channel 
images  or  in  automatic  recognition  maps.  Rabo  images  are  attention-focusing  tools. 

(6)  RAGAL  is  less  expensive  than  statistical  decision  theory  because  SPARC  training  is 
replaced  hv  a  cheaper  training  operation. 

(7)  Because  of  the  relatability  to  labdata,  RAGAL  should  be  capable  of  automatic  train¬ 
ing  in  the  future  — with  lab  data  as  training  sets. 

(8)  Finally,  since  RAGAL  uses  a  much  simpler  logic  than  statistical  decision  theory,  it 
can  be  made  operational  more  easily.  In  fact,  ratio  gating  logic  could  be  done  by  a 
small,  inexpensive  on-board  processor  in  future  aircraft  and  satellites. 

Efforts  are  now  being  made  to  incorporate  the  ERS1S  data  band  into  the  RAGAL  procedure 
and  to  test  the  procedure  again  with  one  or  more  different  geological  data  sets. 


OWME  »L 


WILLOA  RUN  LAHORAIOHitS  THE  UNiVE  WblT  t  Of  MICHIGAN 


Jm 


5 

DISCRIMINATION  OF  COMPOSITIONAL  DIFFERENCES 
BY  THE  RATIO  METHOD 

Aerial  Photography  lias  long  been  an  extremely  useful  tool  for  geologists.  Tie  spectral 
lange  of  photographic  film  sensitivity,  however,  is  limited  to  electromagnetic  spectrum  wave¬ 
lengths  shorter  than  approximately  0.9  pin.  Because  of  this  restricted  spectral  range,  pho¬ 
tography  is  limited  in  two  important  aspects.  First,  the  sun  or  some  other  externa!  source  of 
light  is  always  required  since  only  reflected  radiation  from  a  target  can  be  observed.  Under 
average  day  time  conditions,  with  temperatures  on  the  order  of  300°K,  the  radiation  emanating 
from  targets  at  wavelengths  shorter  than  3  pm  is  primarily  reflective,  and  at  longer  wave¬ 
lengths  is  primarily  emitted  by  the  target  itself.  Secondly,  from  a  geochemical  standpoint, 
photography  is  limited  to  observing  compositional  differences  in  targets  related  to  intra- 

atomic  phenomena.  Electronic  transitions  in  the  outer  shells  of  certain  ions  (such  as  Fe43 
-*2  -*2  -1  ’ 

Fe  ,  Cu  ,  and  OH  )  produce  spectral  reflectance  minima  at  wavelengths  up  to  approximately 

2.3  p m.  However,  there  arc  important  spectral  features  of  silicates  and  other  geological  tar¬ 
gets  (carbonates,  sulfates,  etc.),  related  primarily  to  interatomic  phenomena,  which  occur  at 
wavelengths  longer  than  approximately  2.0  /  m  In  particular,  interatomic  vibrations  of  sili¬ 
cate  rocks  produce  characteristic  spectral  omittance  minima  in  the  8.0-pm  to  14  /.  m  region; 
these  minima  contain  information  concerning  both  the  chemical  composition  and  lattice  struc¬ 
ture  of  crystalline  rocks. 

Over  the  last  decade,  a  new  means  of  collecting  spectral  information  in  imaged  form  has 
been  developed:  the  multispectral  scanner.  Instead  of  detecting  radiation  simultaneously 
from  all  elements  in  the  scene  to  be  imaged  in  only  three  or  fewer  spectral  regions,  as  a 
camera  does,  the  multispectral  scanner  detects  radiation  from  many  narrow  spectral  regions 
for  only  one  point  in  the  scene  at  a  time.  The  final  image  is  produced  through  point -by-point 
scanning.  Finer  spectral  resolution  and  much  greater  spectral  range  (out  to  1-1  /  m  for  ter¬ 
restrial  applications)  are  obtained  with  scanners  than  with  cameras.  Therefore,  the  multi¬ 
spectral  scanner  offers  a  new  and  valuable  source  of  diagnostic  information  about  chemical 
composition  that  can  be  useful  for  geologic  mapping. 

The  first  part  of  this  section  (which  comprises  part  of  Robert  Vincent's  Ph.D.  disserta¬ 
tion  [13  [)  is  limited  to  igneous  silicate  rocks  and  to  the  8 -/nil  to  1 4 -/j  m  spectral  region,  a 
thermal  infrared  atmospheric  "window"  which  contains  some  of  the  most  compositionally 
diagnostic  spectral  information  on  silicates.  Working  with  an  infrared  ratio  method  devised 
earlier  and  a  detailed  knowledge  of  how  detected  radiances  from  a  given  rock  are  related  to 

13  RoberrVincent,  A  Thermal  Infrared  Ratio  Imaging  Method  for  Mapping  Compositional 
Variations  Among  Silicate  Rock  Types,  Unpublished  Ph.D.  dissertation,  The  University  of 
Michigan,  1973. 
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its  chemical  and  mineralogical  compositions,  the  objective  of  this  first  part  of  the  investiga¬ 
tion  is  to  define  a  better  basis  of  discrimination  for  imaging  compositional  differences  among 
silicate  rocks. 

The  second  part  deals  with  experimental  ratioed  imagery  only  in  the  0.4  to  2.5-/im  wave¬ 
length  region,  collected  during  a  1972  flight  over  a  region  in  the  Black  Hills  of  South  Dakota. 
The  rock  types  include  igneous,  metamorphic,  and  sedimentary.  The  objectives  of  this  part 
of  the  investigation  are  to  determine  empirically  which  rocks  in  that  test  area  can  be  dis¬ 
criminated  from  one  another  and  to  develop  a  simpler  method  for  producing  automatic  recog¬ 
nition  maps  of  geologic  targets. 


5.1  DEFINITION  OF  THE  BASIS  FOR  DISCRIMINATION  OF  THE  INFRARED  RATIO  AMONG 

IGNEOUS  SILICATES 

In  last  year's  final  report  under  this  contract  f  14  |  and  in  three  other  publications  |9,  11, 
12 1,  a  two-channel,  infrared-ratio-imaging  techr,..  e  was  described  which  produces  ratio  im¬ 
ages  that  both  enhance  spectral  emittance  variation,  slated  to  compositional  differences  and 
suppress  temperature  variations  across  a  given  scene.  The  infrared  ratio  (called  R.  )  was 

i  ,  C. 

found  to  be  crudely  correlated  (correlation  coefficient  =  0.74)  to  (?,Si02  of  igneous  rocks.  The 
purpose  of  this  part  of  the  investigation  is  to  define  one  or  more  bases  of  discrimination  of 
the  infrared  ratio  method  which  are  better  correlated  than  Si02  with  R}  2> 

In  last  year's  report  [14]  an  experiment  was  described  which  provided  new  laboratory 
data  for  silicate  rocks.  The  spectral  emittance  curves  resulting  from  that  experiment  are 
shown  in  Fig.  18.  The  chemical  composition  of  each  of  the  26  igneous  rocks,  determined 
from  rapid  rock  analysis  by  the  U  S.  Geological  Survey,  is  given  in  Table  5.  From  these 
spectra  and  an  atmospheric  model  [  15 1 ,  the  R  ratio  was  calculated  for  each  rock  for  the 
case  of  the  two-element  infrared  detector  (described  in  last  year's  report  [  1 4  |) flown  at  an 
altitude  of  1  km  above  ground  on  a  clear  drv  day. 

The  equations  for  calculating  Rj  2  are  given  below.  The  total  spectral  radiance  received 
by  the  detector  from  the  direction  of  the  target  can  be  calculated  from  Eq.  (1)  [16 1: 


14.  T.  Wagner,  R.  K.  Vincent,  B.  Drake,  R.  Mitchell,  and  P.  Jackson,  Tunnel-Site  Selec¬ 
tion  by  Remote  Sensing  Techniques,  The  University  of  Michigan,  Technical  Report  No.  10018- 
13-F,  U.S.  Bureau  of  Mines  Contract  H0210041  (ARPA  Order  No.  1579,  Amend.  2,  Program 
Code  IF10),  1972. 

15.  D.  Anding,  R.  Kauth,  and  R.  Turner,  Atmospheric  Effects  on  Infrared  Multispectral 
Sensing  of  Sea  Temperature  from  Space,  Willow  Run  Laboratories,  Report  No.  2676-5-F, 

1970. 

16  W.  L.  Wolfe,  ed.,  Handbook  of  Military  Infrared  Technology,  U.S.  Office  of  Naval  Re¬ 
search,  Department  of  the  Navy,  Washington,  1965,  pp.  9-20. 
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VT)  [f*WT> 


P>Lsky  lrAA  +  LVX 


UL 


bbA^T>  '  LskyA 


■*  L  .  . }  t  . .  L.,  (1) 

skvA1  AA  Va  '  ' 


where 


1^(T)  the  total  spectral  radiance  of  wavelength  a  arriving  at  the  detector 
1  •  the  rock  emittance  taken  from  laboratory  data 
LbbA^T'  the  sPectral  radiance  of  a  blackbody  at  the  same  temperature,  T.  as  the  rock 
LskyA  the  sPect,’al  radiance  from  the  sky  Incident  on  the  diffuse  rock  surface 


rAA  tlie  sPectral  atmospheric  transmittance 

IjVA  "  the  sPertraI  radiance  from  emission  and  scattering  into  the  beam  by  the 
atmosphere  in  the  path  between  rock  and  detector 


For  two  bands  of  approximately  8.2-10.9  pm  and  9.4-12.1  pm  and  an  assumed  rock  tempera¬ 
ture  of  300°K,  Eq.  (D  was  integrated  over  the  respective  spectral  regions,  and  a  ratio  of  ra¬ 
diances  in  the  two  bands  for  each  jf  the  laboratory  rock  samples  was  calculated  from  the 
following  equation: 


Ah 


R 


L1(T)  Aa 

1,2  ~  T^rry =  Ad 


J  fj(A)LA(T)dA 


|  f2(A)LA(T)dA 


Ae 


(2) 


where  Aa  and  Afa  =  limits  of  band  1 
Ac  and  Arf  =  limits  of  band  2 
fj(A)  and  fg(A)  -  the  exact  detector  spectral  responses 
Lj(T)  and  L2(T)  =  the  total  detected  radiances  in  the  two  respective  bands 

This  calculated  R,^  (given  in  Table  6  for  each  rock)  should  agree  with  scanner-measured 
R],2  under  the  altitude  and  atmospheric  conditions  prescribed  above. 

The  question  to  be  answered  is,  What  chemical  and,  or  mineralogical  parameters  of  these 
igneous  rocks  are  best  correlated  with  R>2  ?  Before  the  answers  are  given,  however,  the 
statistical  procedures  used  to  obtain  them ’will  be  explained. 


5  11  STATISTICAL  PROCEDURES 

The  statistical  method  chosen  to  determine  the  best  regression  equations  for  this  investi¬ 
gation  is  a  stepwise  regression  [17|,  which  is  available  in  The  University  of  Michigan’s  MIDAS 


York!7i967,  ^niS”'  Sm'th’  AppHed  RePression  Analysis,  John  Wiley  and  Sons,  Inc..  New 
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TABLE  6.  THERMAL  INFRARED  RATIOS  (R  )  CALCULATED  FROM 

9  ^ 

LABORATORY  SPECTRA 


Sample  Number 

Sample  Name 

Calculated  R 

1 

Granite  (A79) 

1.1761 

2 

Granite  (A119) 

1.1946 

3 

Granite  (A122) 

1.1987 

4 

Rhyolite  (SA-49) 

1.1266 

5 

Rhyolite 

1.1741 

6 

Pyroxene  Syenite  (A83E) 

1.1784 

7 

Trachyte 

1.1442 

8 

Nepheline  Syenite 

1.2094 

9 

Grancdriorite  (A117) 

1.2134 

10 

Granodriorite  (AT27) 

1.1723 

11 

Dacite  (E12A) 

1.1724 

12 

Dacite  (2A1) 

1.1436 

13 

Dacite  (2B5)  W 

1.1574 

Dacite  (2B5)  F 

1.1686 

14 

Diorite  (A129) 

1.1685 

15 

Diorite  (E47C) 

1.1795 

16 

Andesite  (E55A) 

1.1962 

17 

Andesite  (2A3) 

1.1394 

18 

Gabbro  (E80A) 

1.1810 

19 

Basalt  (1A8) 

1.2154 

20 

Basalt  (WI-1-104) 

1.2204 

21 

Basalt  (S95) 

1.2246 

22 

Basalt  (5B4) 

1.2255 

23 

Basalt  (1A3) 

1.2018 

24 

Anorthorite 

1.2458 

25 

Diabase 

1.2204 

26 

Peridotite 

1.3723 
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Sla'lsticai  Laboratory  software  system.  The  stepwise  regression  method  seeks  to  find  the 

best  .near  combination  of  independent  variables.  (X,),  for  prediction  of  the  dependent  variable 
(Y),  through  the  following  steps: 

(1)  Prescribed  significance  levels  are  chosen  for  the  acceptance  and  rejection  of  X 
variables  from  the  regression  equation -for  instance,  a  significance  level  (the* 
higher  the  significance  level  the  more  insignificant  the  variable)  of  0  6  for  accep¬ 
tance  would  mean  that  to  be  selected  for  regression,  an  X  variable  has  to  account 
for  at  least  40^  of  all  the  residual  error  in  Y,  the  predicted  value  of  Y. 

(2)  The  X.  variable  (for  example,  X,)  which  is  most  highly  correlated  with  Y  is  found 

and  a  least  squares  equation  of  the  predicted  value  of  Y  is  calculated,  such  that  Y 

f(X* ), 

(3)  The  remaining  X.  variables  are  searched,  and  the  one  with  the  highest  partial  corre¬ 
late  coefficient  with  Y  (which  will  be  called  X,,)  is  used  to  calculate  a  second  least 
squares  equation,  Y  -  f(Xj.  X^-this  procedure  is  continued  for  the  other  X.. 

(4)  After  each  new  X,  variable  is  selected,  a  partial  F-test  on  the  X.  variables  selected 
in  previous  steps  is  conducted,  and  if  the  significance  level  of  one  of  the  earlier  se¬ 
lected  variables  becomes  higher  than  the  prescribed  significance  level  for  rejection 
the  variable  is  dropped,  and  a  least  squares  equation  is  calculated  without  it -for 
instance,  if  the  addition  of  X3  increases  the  significance  of  X.  above  a  certain  pre¬ 
scribed  level,  the  new  least  squares  equation  would  be  of  the  form  Y  f(X  ,  X  ). 

(5)  The  procedure  terminates  when  there  are  no  other  Xj  variables  with  a  significance 
level  iess  than  the  prescribed  significance  level  of  acceptance. 

This  stepwise  regression  procedure  perform*  primarily  for  equal  .,iB„l„ca„ce 
levels  of  acceptance  and  re|ectio„.  The  procedure  was  performed  once  for  each  of  II, e  fol- 
lowing  significance  levels  of  acceplanee  and  rejection:  0.5,  0.6,  0.7,  0.8,  and  0.99.  The  0  99 
level  accepls  all  of  the  variables  and  rejeels  none  for  almost  any  data  sel. 

Applying  Ihe  slepwise  regression  procedure  for  the  five  prescribed  significance  levels 
resulted  in  from  one  to  five  regression  equations,  each  with  a  different  number  of  variables 
ree  criteria  were  used  to  judge  which  of  these  regression  equalions  was  besl  for  a  given 
se  of  variables.  First,  the  significance  level  of  the  regression  equation  (this  is  different 
from  the  partial  significance  levels  for  X,  variables  mentioned  In  the  above  procedure)  was 
required  to  be  *0.01,  which  mean,  that  there  is  a  99%  or  greater  probability  that  the  result¬ 
ing  regression  equation  was  no,  caused  by  a  se,  of  random  data  points.  Second,  the  multiple 
correlation  coefficient,  fr),  and  standard  error,  (s.e.,-„  the  mean  (V,  value  y-,0r  the 
final  regression  equation  had  ,0  be  improved  over  other  regression  equations  with  fewer 
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variables  to  justify  the  use  of  one  more  variable.  For  n  rock  samples,  r  and  s  e 
by 

fV.-Y,2 

r  i  1 

£v*’2 

j  i 


are  given 


•  V 


Thirdly,  the  regression  equations  had  to  be  physically  reasonable.  The  latter  two  criteria 
were  substantially  subjective. 

Three  of  the  recessions  performed  later  In  this  investigation  Involve  onlv  one  independent 
and  one  dependent  variable.  For  those  eases,  a  simple  correlation  coefficient,  (r  )  will  be 

vnit-vJ.  ’  YV'>  "  UC 


E<v  x,<vi  - 


xv'jvi,  -Jli’CS 

Z<VS  £ 


£v 

>1 


i?on,V!'l7,  o°! V™  ra"E'  rr°'"  '10  '°  ralUCS  of  r  E»-  O)  range  only 


5.1.2  CORRELATION  OF  CALCULATED  RATIOS  WITH  CHEMICAL  PARAMETERS 
The  two  Infrared  channels  for  which  R  have  been  calculated  are  available  on  an  eaisf- 
n, (.infrared  scanner  (the  NASA  M-7  scanner  aboard  the  Environmental  Research  Ins, Rule 
ic  man  s  C-47  aircraft).  The  question  to  be  addressed  in  this  section  is.  With  what  chemical 
parameters  can  a  ratio  of  these  two  channels  be  significantly  correlated?  Because  these 
channels  are  broad  and  overlapping,  the  answers  will  be  pertinent  lo  ratios  for  other,  possibly 
narrower,  channels  in  the  general  vicinity  of  8.2-10.9  pm  and  9.4-12.1  pm.  A  belter  combi¬ 
ne!, on  of  measured  parameters,  given  other  infrared  channels  of  data,  is  determined  in  the 
thesis  mentioned  earlier  ( 13  |  There,  a  linear  combination  of  several  infrared  ratios  is  used 
o  Identify  rock  type  on  fhe  basis  of  a  commonly  used  igneous  rock  classification  chart. 
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In  the  discussion  to  follow,  three  subsets  of  the  data  in  Table  5  will  be  referred  to.  Set  A 
includes  samples  1-25,  excluding  only  26.  Set  B  includes  all  but  four  samples  (1,  3,  5,  and  26) 
These  two  subsets  are  used  only  with  chemical  parameters.  Sample  26,  peridot ite,  was  ex¬ 
cluded  from  consideration  with  the  chemical  parameters  because  it  is  outside  the  range  of 
'<  Si02  for  whlch  a  linear  correlation  between  Si02  and  Rj  2  is  expected  from  previous 
laboratory  data  and  because  its  R^2  ratio  is  sufficiently  distinct  from  the  other  rock  samples 
to  be  unique  Set  C,  used  only  for  the  mineralogical  parameters,  includes  all  samples  except 
for  22  (for  a  reason  given  later). 

Table  7  lists  the  results  of  a  stepwise  regression  performed  on  Set  A  with  oxides  listed 
in  the  order  of  descending  selection.  Similar  results  are  also  shown  for  Set  B.  The  best 
regression  equation  relating  predicted  It j  2  values  (which  uill  be  called  CJ2)  to  12  oxides 
from  Table  5  can  be  expressed  as 


C 


12 


=  A 

o 


12 

y  v, 

i  =  l 


(6) 


where  A^  =  constant  coefficient  (Table  7) 

A.  =  coefficient  for  i-th  oxide  (Table  7) 

O.  =  weight  y  of  i-th  oxide  (e.g.,  ^Si02,  ((  Fe^,  etc.) 

For  Set  A,  the  correlation  coefficient,  (r),  between  the  chemical  parameter  C12  from  Eq.  (6) 
and  experimental  R]2  values  from  Table  2  is  0.9131,  with  an  approximate  0.018  standard 
error  in  R^.  Part  of  the  reasoning  which  led  to  the  use  of  all  12  oxides  was  based  on  the 
fact  that  a  geologist  who  has  a  rapid  rock  analysis  performed  on  his  samples  can  have  all 
12  of  these  oxides  measured  as  cheaply  as  it  would  cost  him  to  have  only  nine  or  ten  mea¬ 
sured. 


The  best  single  chemical  parameter  is«(SiC>2,  which  is  plotted  versus  Rj  2  in  Fig.  19 
for  all  26  rock  samples.  The  simple  correlation  coefficient  between  f'(SiC>2  and  R  for  Set 
A  is  -0.74,  with  a  standard  error  in  Si02  of  6.39,  where  Rj  g  is  taken  to  be  the  independent 
variable.  From  Fig.  19  it  is  obvious  that  samples  1,  3,  and  5  are  responsible  for  much  of 
the  error.  Hence,  a  regression  between  f?Si02  and  Rj  2  was  performed  for  Set  B  (excluding 
samples  1,  3,  and  5  as  well  as  26),  with  the  resulting  simple  correlation  coefficient  of  -0.87 
and  standard  error  in  f;(Si02  of  4.06,  where  Rj  2  is  once  again  the  independent  variable. 
Samples  1  and  3  are  coarse  crystalline  granites.  Part  of  the  reason  they  are  atypical  may 
be  their  large  grain  size  in  relation  to  the  spot-size  used  for  the  infrared  laboratory  measure 
men's  (as  discussed  in  the  previous  section).  Sample  5  is  a  rhyolitic  welded  tuff.  It  cannot 
be  conclusively  determined  from  this  data  set  whether  these  samples  are  atypical  for  natural 
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TABLE  7.  REGRESSION  OF  THE  C^  (SET  A)  and  C  (SET  B) 
CHEMICAL  PARAMETERS  WITH  R 

1 1  ^ 


Oxide 

Same 

No. 

Cj  2 (Set  A) 
Coef  f lc ients 

Cj 2 (Set  B)  | 

Coef  fic ients 

SiO,, 

1 

-.000153 

-.002149 

F«2»3 

2 

-.020216 

-.017907 

H2° 

3 

-.000126 

.004483 

Na20 

4 

.009202 

.009342 

CaO 

5 

.007289 

.006034 

1  2° 

6 

.005383 

.005443 

MgO 

7 

.008184 

.004965 

no2 

8 

-.010045 

-.019345 

FeO 

9 

-.002312 

-.001141 

MnO 

10 

.050832 

.030810 

P2°5 

11 

.015214 

.036184 

A12°3 

12 

.000740 

.000311 

Constant 

0 

1.111000 

1.242200 

Correlation 

(r) 

.91 

.95 

Standard 

Error (s. e. ) 

.018 

.016 
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reasons  or  for  reasons  related  to  laboratory  measurements.  Nevertheless,  the  relatively 
high  correlation  between  SiC>2  and  Rj  2  t°r  the  remaining  22  silicate  rocks  (excluding 
peridotite)  is  significant. 

The  second  chemical  parameter  chosen  by  the  stepwise  regression  was'.  Fe2Og,  such 
that  for  Set  A  the  equation  for  the  predicted  value  of  Rj  2  on  the  basis  of',  SiC>2  and  '>  Fc2®3 
is  given  by 

R]  1.3984  -  0.003069(',Si02)  -  0.016486(f;  Fc203)  1  3984  -  0.003  069  {'y  |Si02  + 

5.4(Fe203)|[ 

Tliis  indicates  that  the  best  two-oxide  parameter  combination  is  approximately  ,r(  [SiC>2  + 

5  4 ( Fe2Og )  | .  This  parameter  is  plotted  versus  Rj  2  in  Fig.  20.  The  simple  correlation  of 
%[SiO„  +  5.4(Fe„0„)|  versus  R.  „  is  -0.82,  with  a  standard  error  in  the  two-oxide  parameter 
of  4.63,  where  Rj  2  is  the  Independent  variable.  Similar  results  were  obtained  for  Set  B,  ex¬ 
cept  that  the  coefficient  of  rr(  Fe2®3  produced  the  highest  correlation  was  3.7.  The  new 

parameter  is  plotted  in  Fig.  21  for  members  of  Set  B  only.  The  simple  correlation  for  Set  B 
is  -0.92,  with  a  standard  error  in  [Si02  +  3.7(Fc203)|  of  2.94,  where  Rj  2  is  again  the  inde¬ 
pendent  variable. 

In  Table  7,  r  ’,  A1,,03  is  the  last  chemical  parameter  to  be  selected  by  the  forward  regres¬ 
sion  procedure.  However,  if Fe203  is  omitted  from  the  list  and  a  regression  is  performed 
on  the  remaining  11  chemical  parameters,  ',  AlgOg  is  the  second  parameter  selected.  The  co¬ 
efficients  of  'VSi02  and  AlgOg  are  almost  equal  but  of  opposite  sign  for  this  case.  Hence, 
the  parameter  °}  (SiC>2  -  AlgOg),  which  is  plotted  versus  Rj  2  for  both  Set  A  and  Set  B  in  Fig. 
22,  is  a  better  parameter  than^SiO,,  and  only  slightly  worse  than  'J  [SiC>2  +  3.7(Fe203)|  for 
Set  B.  The  simple  correlation  for  Set  A  is  -0.76,  with  a  standard  error  of  7.43  in  (SiC>2  - 
AlgOg):  for  Set  B  it  is  -0.90,  with  a  standard  error  of  4.24  in  "J  (Si02  -  AlgOg). 

The  fact  that  ",  Fe203  should  be  such  an  important  oxide  for  relating  the  chemical  and 
infrared  emittance  properties  of  silicate  rocks  is  difficult  to  understand,  since  Fe2C>3  itself 
exhibits  few,  if  any,  spectral  features  in  the  8-  to  12-pm  wavelength  region.  Ferric  oxide 
(Fe2C>3)  is  a  late  differentiate  of  magma,  as  well  as  a  product  of  the  weathering  of  mafic 
minerals.  However,  since  the  coefficients  of  Fe203  and';;,Si02  have  the  same  sign,  the 
effect  of  more  Fe,,03  is  the  same  as  the  effect  of  higher  ^Si02  (i.e.,  an  increase  in  Rj  2  — 
in  the  direction  of  more  acidic  rocks).  This  may  indicate  that  the  late  differentiation  of  Fe203 
is  more  important  with  regard  to  infrared  spectral  features  than  the  fact  that  ferric  oxide  is 
also  a  weathering  product  of  mafic  minerals.  The  H20,  also  associated  with  late  differentia¬ 
tion,  Is  also  negatively  correlated  with  R.  _.  There  is  a  possible  correlation  between  %  Fe  O 
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FIGURE  21.  CORRELATION  OF  <£Si02  -  3.7(%Fe203)  WITH  Rj  2  FOR  SET  B 


FIGt'RF  22  CORRELATION  OF'y  SiC^  -  <]}  A1203  WITH  R,  ,  FOR  SET  A  AND  SET  B 
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and  the  texture  of  acidic  rocks,  as  evidenced  by  the  fact  that  the  average  Ft'g Og  >s  approxi¬ 
mately  twice  as  high  for  the  acidic  fine-  and  medium-grained  rocks  (samples  4,  5,  7.  and  17) 
as  it  is  for  the  coarse-grained  acidic  rocks  (samples  1,  2,  3,  and  10)  —according  to  the  new 
rock  classifications  in  Table  9  of  the  next  section.  The  intermediate  and  basic  rock  samples 
show  about  the  same  Fe^O^  for  coarse-  and  fine-grained  rocks.  Hence,  the  importance  of 
Fp2®3  *n  ^12  could  possibly  be  an  indication  of  some  textural  effect  on  the  Rj  2  ratio  for 
acidic  rocks.  The  correlation  of  greater  ' ,  Fo2C>3  smaller  grain  size  in  acidic  rocks 
could  be  related  either  to  the  near-surface  origin  of  the  magma  (where  sedimentary  Pe-jO^ 
could  enrich  the  magma)  or  to  disseminated  thin  coatings  formed  preferentially  on  small 
crystals.  The  fact  that  the  intermediate  and  basic  rock  samples  showed  no  such  correlation 
is  evidence  contrary  to  the  latter  hypothesis, 

The  relatively  high  correlation  of  f,'(Si02  -  A 1  g )  with  2  is  not  so  surprising,  since 

the  infrared  emittance  minima  (reststrahlen  bands)  of  feldspars  occur  at  longer  wavelengths 

^3  *4 

than  quartz.  The  substitution  of  A1  for  Si  in  SiO^  tetrahedra,  which  reduces  the  number 
of  Si-0  bonds  and  increases  the  Al-O  bonds,  is  common.  From  our  earlier  discussion 
concerning  the  lower  wavelength  position  of  the  Si -O  reststrahlen  bands  relative  to  those 
caused  by  Al-0  bonds,  it  is  understandable  that  increased  A1203  would  cause  the  emittance 
minima  to  occur  at  longer  wavelengths,  thereby  increasing  Rj  2-  Although  the  selection  of 
Fp2®3  *n  t,ie  fonvarcl  regression  renders  the  r,  AlgOg  insignificant,  ferric  oxide  is  chosen 
second  to  (”( (Si02  -  A1,,03)  in  forward  regression  if (Si 02  -  A1203)  is  treated  as  one  param¬ 
eter.  Hence,  the  utilization  of  Fe203  suppresses  the  importance  of  A,2°3>  but  t,ie  con" 
verse  is  not  true. 

Of  these  chemical  parameters,  Cj2  is  the  best  for  interpreting  R^  2  scanner  ratio  images, 
if  the  geologist  has  a  collection  of  rocks  for  which  the  rapid  rock  analyses  are  known  and  if 
he  has  no  normative  calculation  program  to  convert  the  oxides  into  mineralogical  variables, 
or  no  modal  analyses  of  the  rocks.  The  CJ2  parameter  can  be  calculated  for  his  collection 
from  Eq.  and  Table  7,  and  the  resulting  Cj2  values  can  be  compared  with  Rj  2  from 
scanner  data.  If  either  a  normative  calculation  computer  program  or  modal  analysis  results 
are  available,  however,  the  mineralogical  parameters  described  in  the  next  section  are 
superior  to  CJ2  for  that  purpose. 

The  chemical  parameters  a}  (SiC>2  -  AlgOg),  9,  [SiC>2  +  5.4(Fe203)|,  and  C12  for  igneous 
rocks  may  be  more  relevant  than  mineralogical  parameters  for  classifying  soils  by  the  R. 

J  ,  C. 

ratio  imaging  method,  because  igneous  rocks  and  soils  are  mineralogically  different  but 
contain  similar  oxides.  The  (Si02  -  AlgOg )  parameter  would  seem  to  have  the  best  chance 

of  these  three  for  being  applicable  to  both  soils  and  igneous  rocks,  because  the  replacement 

+3  +4 

of  A1  for  Si  is  a  common  effect  of  weathering. 
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!>.l  3  CORRELATION  OF  CALCULATED  RATIOS  WITH  MINFRALOGICAL.  PARAMETERS 

Instead  of  stopping  with  the  C^  chemical  parameter,  for  two  reasons  it  is  also  desirable 
to  relate  Rj  2  ratios  with  the  volume  percentage  of  minerals  present  in  the  rock  samples. 
First,  since  the  infrared  spectral  features  in  the  8-  to  14-pm  wavelength  region  arise  from 
intramolecular  energy  level  transitions,  mineral  parameters  should  describe  the  molecules 
present  better  than  chemical  parameters  can.  For  instance,  <}  Si02  does  not  represent  just 
the  free  quartz  in  a  rock;  it  includes  silicon  and  oxygen  from  other  minerals  that  need  not 
even  have  teetosilicate  structures.  Second,  quantitatively  accurate  oxide  percentages  must 
be  obtained  through  rapid  rock  analysis,  which  is  moderately  expensive  and  slow  for  geologists 
who  do  not  have  their  own  laboratory  facilities  for  wet  chemical  analysis.  It  is  desirable, 
therefore,  from  the  standpoints  of  cost  and  convenience  alone  to  relate  the  R  ratios  to 
relative  amounts  of  various  common  rock -forming  minerals,  most  of  which  are  identifiable  by 
visual  inspection  and  microscopy.  Volume  instead  of  weight  percentages  are  desirable  be¬ 
cause  the  former  are  more  relatable  to  area  percentages  in  a  thin  section  and  a  rock  surface 
imaged  by  a  multispectral  scanner. 

The  mineralogical  volume  percentages  for  this  experiment  were  calculated  from  a 
C  ?  P  w-  norin  computer  program  of  Rodgers,  Cochrane,  and  LeCouteur  f  18  J,  on  the  basis 
of  the  oxide  percentages  obtained  from  raptd  rock  analysis  (see  Table  5).  The  C.I.P.W.  rock 
classification  method  was  created  by  Cross,  Iddings,  Pirsson,  and  Washington  [19 1.  It 
classified  rocks  accordtng  to  the  proportions  of  certain  normative  mtnerals,  whtch  are  listed 
in  Table  8.  Thes->  are  not  all  the  same  as  the  modal  minerals  teat  actually  occur  in  rocks, 
though  some  are  ,he  same  for  both  cases.  The  norm  is  the  calculated  mineral  composition 
of  a  given  roc!.,  as  opposed  to  the  actual  mineralogtcal  composition  or  mode.  Arbitrary 
divisions  [20 1  are  made  from  the  proportions  of  the  normative  minerals.  The  inputs  to  the 
C.I.P.W.  classification  scheme  are  chemical  oxides  of  a  given  rock  such  as  those  shown  in 
Table  5,  and  the  output  is  the  calculated  mineral  composition  of  the  rock,  in  terms  of  the 
normative  minerals. 

Table  9  gives  the  mineralogical  parameters  resulting  from  the  Rodgers,  Cochrane,  and 
LeCouteur  computer  program  (modified  to  yield  volume  percentages)  for  Set  C.  Sample  22, 

18.  K.  A.  Rodgers,  R.  H.  A.  Cochrane,  and  P.  C.  LeCouteur,  Fortran  II  and  Fortran  IV 
Programs  for  Petrochemical  Calculations,  Mineralogical  Magazine,  37,  1970,  pp.  952-953. 

19.  W.  Cross,  J.  P.  Iddings,  L.  V.  Pirsson,  and  H.  S.  Washington,  A  Quantitative  Chemico- 
Mineralogical  Classification  and  Nomenclature  of  Igneous  Rocks,  J.  of  Geology,  10,  1902,  pp. 
555-690 

20.  A.  Johannsen,  A  Descripttve  Petrography  of  the  Igneous  Rocks,  Vol.  I,  The  University 
of  Chicago  Press,  Chicago,  1931,  pp.  85-88. 
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TABLE  8.  NORMATIVE  MINERALS 


Abbreviation 


Ounrtz  OTZ 

Or thoclase  ORTHO 

Aibite  ALB 

Anorthite  ANOR 

Nepheline  NEPH 

Clinopvroxene  (Calcium-rich)  CP (WO) 

Clinopyroxene  (Magnesium-rich)  CP(EN) 

Clinopyroxene  (Iron-rich)  CP(FS) 

Orthopyroxene  (Magnesium-rich)  OP(EN) 

Orthopyroxene  (Iron-rich)  OP(FS) 

Wollastonite  W0L 

Olivine  (Forsterite)  OL(FO) 

Olivine  (Fayalite)  OL (FA) 

Hematite  HEM 

Magnetite  MG 

Ilmenite 

Corundum  GqE 

Apatite  ApA 

Caicite  CAL 

Water  WAT 
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TABLE  9.  COMPUTED  "NORMS"  OF  SET  C  ROCK  SAMPLES 
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a  basalt,  was  omitted  because  its  calculated  norm  included  a  normative  mineral  which  is  not 
found  in  nature,  a.'.d  since  there  were  enough  other  basalts  to  represent  that  rock  type  ade¬ 
quately. 

There  were  two  major  reasons  why  the  norms  were  used  instead  of  modes  as  mineralogical 
parameters.  First,  it  is  difficult  to  perform  quantitatively  accurate  modal  analysis  on  aphanitic 
ground  masses,  which  are  present  in  some  of  the  samples.  The  resulting  modes  can,  therefore, 
be  more  descriptive  of  the  phenocrysts  than  of  the  ground  mass.  The  normative  calculation,  by 
contrast,  is  influenced  equally  by  phenocrysts  and  ground  mass,  because  the  chemical  oxides 
on  which  the  calculation  is  based  are  bulk  oxides  obtained  from  wet  chemical  analysis  of 
crushed  samples.  Since  an  infrared  scanner  or  spectrometer  looks  at  both  phenocrysts  and 
ground  mass,  the  bulk  properties  should  be  better  correlated  with  the  infrared  spectral  fea¬ 
tures.  Secondly,  since  the  normative  calculations  of  Table  9  are  derived  solely  from  the  ox¬ 
ides  of  Table  5,  the  correlations  of  mineralogical  and  chemical  parameters  with  R.  „  can  be 

1 

directly  compared,  with  little  probability  of  unequal  treatment  because  of  mineral  mis-identi- 
fication  or  mis-estimate  of  volume  percentage. 

The  normative  computer  program  was  found  to  agree  well  with  a  separate  normative  pro¬ 
gram  of  N.  G.  Ware  [21 1  for  calculated  norms  of  several  basic  rock  samples  [22  | .  The  norm 
of  a  rock  may  not  closely  agree  with  its  mode,  especially  for  pyroxene-rich  rocks  [ 22 1 .  How¬ 
ever,  the  norms  and  modes  usually  agree  within  a  few  percent  for  most  igneous  rocks,  within 
the  context  of  the  summations  implicit  to  the  normative  minerals  (e.g.,  orthoclase  encompass¬ 
ing  all  potassic  feldspars). 

Stepwise  regressions  were  performed  between  Rj  2  and  the  computer  norms  of  Table  9, 
which  represent  all  of  the  samples  in  Set  C.  The  best  regression  equation  coefficients  are 
given  in  Table  10.  The  coefficients  in  Table  10  can  be  used  to  establish  a  mineralogical  index 
defined  by  the  following  equation: 

16 

M.6-Bo*I>iMi  «« 

i  =  l 

where  Mjg  =  mineralogical  index  for  16  mineral  parameters 
Bq  =  constant  coefficient  (Table  10) 

B.  -  coefficient  for  i-th  mineral  parameter  (Table  10) 

M.  =  volume  ^  of  i-th  normative  mineral 

21.  N.  G.  Ware,  Department  of  Geophysics,  Australian  National  University,  Canberra, 

Private  Communication,  1973. 

22.  E.  Essene,  Department  of  Geology  and  Mineralogy,  The  University  of  Michigan,  Ann 
Arbor,  Private  Communication,  1973. 


71 


1 


*  OHMI  BLY  Will  OW  HUN  t  A  IIOMA  I  OHlt  b  1MI  UMVf  HSDr  Ol 


TABLE  10.  REGRESSION  OF  THE  M. ,  MINERALOGICAL  INDEX 

lo 

WITH  Rx  2  FOR  SET  C 


Mineral 

Parameter  No. 

Coefficient 

Olivine  (Forsterite) 

1 

-.002417 

Wollastonite 

2 

. 350000 

Water 

3 

-.013110 

Hematite 

4 

-.112930 

Albite 

5 

-.002834 

Quartz 

6 

-.001032 

Olivine  (Fayalite) 

7 

.022932 

Ilmenite 

8 

-.090608 

Apatite 

9 

.120840 

Orthopyroxene 

(Magnesium  Rich) 

10 

.005732 

Nepheline 

11 

-.006352 

Corundum 

12 

.008251 

Clinopyroxene 
(Iron  Rich) 

13 

.189930 

Clinopyroxene 
(Calcium  Rich) 

14 

-.129250 

Clinopyroxene 
(Magnesium  Rich) 

15 

.128740 

Orthopyroxene 
(  Iron  Rich) 

16 

-.011553 

Constant 

0 

1.3606 

Correlation^ ) 

.99 

Standard  Error 
(s.e. ) 

.009 
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The  Mjg  niineralogical  index  is  actually  the  predicted  value  of  R  based  on  the  regression 

represented  by  Table  10.  The  multiple  correlation  coefficient  between  and  measured  values 

lb 

Of  Rj  2  from  Table  7  is  0  99,  with  a  standard  error  of  0.009  in  Rj  2  —  if  Mjg  is  used  to  predict 
the  measured  Rj  2 •  A  plot  of  Mjg  versus  Rj  2  is  shown  in  Fig,  23. 

Up  to  this  point,  the  names  associated  with  the  rock  samples  have  been  those  assigned  by 
their  respective  collectors,*  who  named  them  on  the  basis  of  either  hand  specimen  inspection 
or  modal  analysis,  and  probably  according  to  slightly  different  schemes.  In  an  attempt  to  place 
these  rock  classifications  on  a  more  uniform  basis,  each  rock  was  reclassified  according  to 
Travis  s  classification  of  igneous  rocks  [23  |,  on  the  basis  of  the  volume  percentage  of  norma¬ 
tive  minerals  given  in  Table  9.  Whenever  the  rock  was  a  borderline  case  m  this  classification 
scheme,  the  collector’s  name  for  the  specimen  was  retained.  Table  11  shows  the  resulting  new 
sample  names,  along  with  the  old  names,  arranged  in  order  of  increasing  Rj  2.  Sample  22,  not 
a  member  of  Set  C,  is  also  given,  with  unchanged  classification.  In  the  ncxt-to-last  column, 
a  group  name  is  assigned  to  each  sample,  according  to  Travis's  nomenclature.  The  group 
names  are  explained  in  Table  12. 

The  most  important  point  to  be  made  from  Table  11  is  that  (as  shown  by  the  last  two  columns 
of  Table  9)  the  Rj  2  ratio  is  not  well  correlated  with  the  rock  groupings  of  Travis's  nomencla¬ 
ture.  On  this  basis,  only  the  I  and  ,T  groups  seem  to  be  identifiable  in  the  presence  of  other 
silicates.  This  brings  out  an  important  point:  the  minerals  on  which  geologists  have  chosen  to 
place  greatest  importance  in  their  traditional  rock  classification  schemes  do  not  have  the  same 
relative  importance  with  regard  to  their  effect  on  the  infrared  properties  of  rocks.  For  instance, 
Travis's  classification  places  greatest  importance  on  quartz,  potassic  feldspar,  feldspathoids, 
pyroxene,  and  olivine,  with  secondary  importance  on  the  relative  amounts  of  sodic  and  calcic 
plagioclase.  However,  when  a  stepwise  regression  was  performed  to  find  the  eight  minerals 
most  significant  for  predicting  the  infrared  ratio,  the  results  were  forsterite,  anorthite,  wol- 
lastonite,  water,  magnetite,  hematite,  albite,  and  quartz  —  in  that  order.  Since  water,  magne¬ 
tite,  and  hematite  have  almost  no  infrared  spectral  features,  they  must  influence  the  infrared 
ratio  in  some  indirect  manner.  Quartz  and  forsterite  have  reststrahlen  bands  at  the  shortest 
(8.1-10.1  pm)  and  longest  (10.2-12.2  pm)  wavelengths  of  all  the  silicates.  Albite  and  anorthite 
reststrahlen  bands  fall  at  9.6-10,1  pm  and  8.7-10,9  pm,  respectively.  Wollastonite  (CaSiO^), 
a  pvroxenoid  with  a  1:3  silicon  to  oxygen  ratio,  is  expected  to  have  reststrahlen  bands  just 


•These  rocks  were  collected  by  several  geologists  and  loaned  to  Dr.  Lawrence  Rowan  of 
the  U.S.  Geological  Survey,  who  in  turn  loaned  them  to  Robert  Vincent. 

23.  R.  B.  Travis,  Classification  of  Rocks,  Quarterly  of  the  Colorado  School  of  Mines  50 
1955,  p  12  ’  ’ 
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1.40 


FIGURE  23.  MINERALOGICAL  PARAMETER  MJ6  VERSUS  Rj  2  FOR  SET  C 
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1ALL^  J1;  iRECLASSIFICATIOK  of  rock  samples  according  to 
c.i.p.w.  calculations  and  iravis's  nomenclature 


Sample 

No. 

Sample 

Old  Name 

Sample 

New  Name 

Assign 

Group 

4 

Rhyolite(SA-49) 

Quartz  Latite 

D-3 

17 

Rhyolite (2A3) 

Rhyolite 

A- 3 

12 

Dacite(2Al) 

Dacite 

G-3 

7 

Trachyte 

Quartz  Latite 

D-3 

13 

Dacite(2B5)W 

Dacite 

G-3 

14 

Diorite (A129) 

Granodiori te 

G-l 

10 

Granadiorite  (A127] 

1  r  lartz  Monzonite 

D-l 

11 

Dacite (E12A) 

uartz  Latite  Porphyry 

D-2 

5 

Rhyolite  (Welded 
Tuff) 

Quartz  Latite 

D-3 

1 

Granite(A79) 

Quartz  Monzonite 

D-l 

6 

Pyroxene  Syenite 
(E83A) 

Pyroxene  Syenite 

B-l 

15 

Diorite (E47C) 

Diorite 

H-l 

18 

Gabbro(E80A) 

Gabbro 

1-1 

2 

Granite(A119) 

Quartz  Monzonite 

D-l 

16 

Andes ite(E55A) 

Dacite 

G-3 

3 

Granite(A122) 

Quartz  Monzonite 

D-l 

23 

Basalt (1A3) 

Basalt 

1-3 

8 

Nepheline  Syenite 

Nepheline  Monzonite 

F-l 

9 

Granodiorite 
(All  7) 

Diorite 

H-l 

19 

Basalt(lA8) 

Basalt 

1-3 

25 

Diabase 

Diabase 

1-1 

20 

Basalt (WI-1-104) 

Basalt 

1-3 

21 

Basalt (S95) 

Basalt 

1-3 

22 

Ba6alt(5B4) 

Basalt 

1-3 

24 

Anorthosite 

Anorthosite 

1-1 

26 

Peridot ite 

Peridotite 

J-l 

Ratio 

1.1266 

1.1394 

1.1436 

1.1442 

1.1574 

1.1685 

1.1723 

1.1724 

1.174] 

1.1761 

1.1784 

1.1795 

1.1810 

1.1946 

1.1962 

1.1987 

1.2018 

1.2094 

1.2134 

1.2154 

1.2204 

1.2204 

1.2246 

1.2255 

1.2458 

1.3723 
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TABLE  12.  ASSIGNMENT  OF  GROUP  CLASSIFICATIONS  ACCORDING  TO 

TRAVIS'S  NOMENCLATURE 


Symbo1  Interpretation 


A 

B 

C 

D 

E 

F 

G 

H 

I 


J 

1 

2 

3 


Potassic  Feldspar  >2/3  Total  Feldspar,  Quartz  >10% 

P°£aiHiC  >2/3  T°tal  Feldspar,  Quartz  <10%, 

Feldspathoid  <10% 

Potassic  Feldspar  >2/3  Total  Feldspar, 

Feldspathoid  >10% 


Potassic  Feldspar  1/3  to  2/3  Total 

Potassic  Feldspar  1/3  to  2/3  Total 
Feldspathoid  <10% 


Feldspar , 
Feldspar , 


Quartz  >10% 
Quartz  <10%, 


Potassic  Feldspar  1/3  to  2/3  Total 
Feldspathoid  >10% 


Feldspar , 


Plagioclase  >2/3  Total  Feldspar, 
Total  Feldspar,  Quartz  >10% 


Potassic  Feldspar  >10% 


Flagioclase  >2/3  Total  Feldspar,  Potassic  Feldspar  <107 
Totai  Feldspar,  Quartz  <10%,  Feldspathoid  <10%, 

Sodic  >  Calcic  Plagioclase 


Plagioclase  >2/3  Total  Feldspar 
Total  Feldspar,  Quartz  <10%, 
Calcic  >  Sodic  Plagioclase 


,  Potassic  Feldspar 
Feldspathoid  <10%, 


<10% 


Little  or  No  Feldspar,  Chiefly  Pyroxene  and/or  Olivene 
Phaneritic  Texture 


Porphyritic  Texture 
Aphanitic  Texture 
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short  of  those  for  anorthite,  since  both  contain  Ca-O  bonds.  In  the  stepwise  regression  for 
Mfg-  anorthite  and  magnetite  were  eventually  replaced  by  other  minerals 

The  fact  remains  that  R  „  discriminates  among  silicates  on  the  basis  of  the  M,,  miner- 

i  *  16 

alogical  parameter,  which  is  different  from  the  basis  of  traditional  rock  classification  schemes. 
Table  13  shows  the  M16  values,  in  ascending  order,  for  the  rock  samples  of  Set  C.  Also  shown 
in  Table  13  is  a  suggested  classification  scheme  based  on  values  of  MJ(,.  Category  III  contains 
mostlj  basalts  and  diabase,  Category  IV  is  anorthosite  and  possibly  some  other  very  basic 
rocks,  and  Category  Vis  strictly  ultrabasic.  Categories  I  and  II  are  not  well  correlated  with 
rock  groups  based  on  Travis's  nomenclature,  except  that  group  I  contains  five  of  the  seven 
fine-grained  volcanic  samples.  This  could  be  indicative  of  a  textural  effect  on  Rj  2  for  acidic 
rocks.  However,  this  effect  could  be  indirectly  caused  by  composition — for  example,  the  D-3 
and  G-3  samples  in  group  I  all  have  higher  %  water  than  the  D-l  and  G-l  samples,  respectively. 

The  value  of  Eq.  (8)  is  that  if  a  geologist  knows  the  norms  of  a  collection  of  his  rocks,  he 
can  calculate  MJ6  (which  is  really  the  best  predicted  value  of  Rj  2)  without  having  to  make  in¬ 
frared  laboratory  measurements  on  the  rocks.  If  the  detector  which  measures  Rj  2  is  flown 
at  1  km  above  the  ground  over  his  test  site,  the  Rj  2  values  measured  by  the  scanner  can  then 
be  used  to  help  him  decide  which  of  the  rocks  in  his  collection  a  given  outcrop  (seen  on  the 
Rl,2  ratio  image)  could  or  could  not  be.  Given  an  Rj  2  value,  the  problem  of  determining  the 
standard  deviations  in  is  a  difficult  statistical  one  involving  inverse  regression  methods 
for  several  independent  parameters  (the  mineral  volume  percentages  in  this  case).  This 
problem  has  not  been  solved,  and  it  was  not  the  subject  of  any  published  statistical  papers 
between  1961-69  [24 j.  To  estimate  crudely  this  standard  deviation  in  MJ6  (or  predicted 
value  of  R1j2).  Biven  the  value  of  Rj  2,  one  can  compare  the  MJ6  values  from  Table  13  with 
the  R^2  values  of  Table  11.  The  largest  disagreement  between  MJ6  and  Rj  2  is  0.011.  There¬ 
fore,  the  standard  deviation  from  predicting  MJ6  from  measured  Rj  2  can  be  estimated  as 
approximately  Oj  =  0.01.  The  estimated  standard  deviation  from  an  error  in  the  measurement 
of  R^2  is  approximately  a 2  =  0.007.  The  combined  standard  deviation,  assuming  that  additive 
errors  can  be  expressed  as  in  Wilson  [25  |,  is 

/  n  o 

o  =  1°!  +  °2  =  0  012  (9) 


24.  B.  L.  Joiner,  N.  F.  Laubscher,  E.  S.  Brown,  and  B.  Levy,  AnAuthorand  Permuted 
Title  Index  to  Selected  Statistical  Journals  (1961-1969),  National  Bureau  of  Standards  Publica¬ 
tion  321,  U.S.  Department  of  Commerce,  Washington,  1970. 

25.  E.  B.  Wilson,  Jr.,  An  Introduction  to  Scientific  Research,  McGraw-Hill  Book  Co  Inc 

New  York,  1952,  p.  273.  ’’  '* 


77 


^ERJM _ 

fO»«t»L»  ftILlUW  HUNlABriMlollits  THl  UNvlHSUYOt  MK_HK,»N 


1  ABLE  13.  A  CLASSIFICATION  BASED  ON  THE  MINERALOGICAL  INDEX 


Sample 

No. 

Old  Sample 

Name 

New  Sample 

Name 

M 

16 

Travis's  Groups  Groups 

V 

4 

Rhyolite 

Quartz  Latite 

D-3 

1.1255 

13 

Dac i te 

Dacite 

G-3 

1.1412 

17 

Rhvolite 

Rhyolite 

A- 3  I  (<L.160) 

1.1423 

12 

Da.'ite 

Dacite 

G-3 

1.1555 

7 

Trachyte 

Quartz  Latite 

D-3 

1.1578 

10 

Granodiorite 

Quartz  Monzonite 

D-l 

1.1691 

14 

Diorite 

Granodiorite 

G-l 

1.1702 

5 

Rhyolite 

Quartz  Latite 

D-3 

1.1708 

11 

Dacite 

Quartz  Latite  Porphyry 

D-2 

1.1726 

1 

Granite 

Quartz  Monzonite 

D-l 

1.1740 

6 

Pyroxene  Syenite 

Pyroxene  Syenite 

B-l  11(1.160-1.200) 

1.1741 

15 

Diorite 

Diorite 

H-l 

1.1793 

18 

Gabbro 

Gabbro 

1-1 

1.1823 

2 

Granite 

Quartz  Monzonite 

D-l 

1.1931 

16 

Andesite 

Dacite 

G-3 

1.1996 

23 

Basalt 

Basalt 

1-3 

1.2033 

8 

9 

Nepheline  Syenite 

Granodiorite 

Nepheline  Monzonite 

Diorite 

F-l 

H-l 

1.2094 

1.2115 

19 

Basalt 

Basalt 

1-3111(1.200-1.230) 

1.2146 

20 

Basalt 

Basalt 

1-3 

1.2181 

25 

Diabase 

Diabase 

1-1 

1.2188 

21 

Basalt 

Basalt 

1-3 

1.2271 

24 

Anorthosite 

Anorthosite 

1-1  IV(1. 230-1. 270) 

1.2446 

26 

Peridot ite 

J-l  V (>1.270) 

1.3724 
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Twice  this  standard  deviation  is  0.024,  which  represents  the  minimum  separation  in  M^g  value 
that  two  rocks  must  have  in  order  to  be  discriminated  with  a  68%  probability  of  success  by  a 
measurement  of  R^  ^ 

Although  the  index  is  the  best  single  parameter  correlated  with  Rj  2  (until  the  sta¬ 
tistical  base  of  rock  samples  is  extended  beyond  the  25  samples  of  Set  C),  the  geologist  may 
find  that  a  mineralogical  index  that  is  more  simply  determined  is  more  useful  with  rock 
samples  for  which  he  has  no  rapid  rock  analyses.  In  that  case,  it  would  be  helpful  to  have  a 
mineralogical  parameter  based  only  on  minerals  recognizable  with  a  hand  lens  or  polarizing 
microscope.  Table  14  shows  the  results  for  Set  C  of  confining  the  mineralogical  parameters 
to  quartz,  anorthite,  nepheline,  olivine  (forsterite  plus  fayalite),  magnetite,  and  orthoclase. 
Because  of  the  nature  of  C.I.P.W.  calculations,  anorthite  represents  all  calcic  plagioclase  and 
orthoclase  represents  all  potassic  feldspars.  Olivine  for  Table  14  was  taken  to  be  sum  of 
forsterite  and  fayalite.  The  resulting  simplified  mineralogical  index  is  as  follows  (where  % 
refers  to  volume  % ): 

Mg  =  1.2081  -  0.0008(%  Quartz)  +  0.0007(%  Calcic  Plagioclase)  +  0.0027(%  Nepheline) 

+  0.0023(%  Olivine)  -  0.0169(%,Magnetite)  -  0.0005(%  Potassic  Feldspar)  (10) 

The  correlation  between  Mg  and  Rj  2  is  0.94:  the  standard  error  in  Rj  2  when  the  latter  is 
the  independent  variable  is  0.019.  The  estimated  combined  standard  deviation  in  Mg  when 
Rj  2  is  known  (calculated  in  the  same  manner  as  that  for  Mjg  in  Eq.  (9)  is  0.030.  Though 
Mg  is  more  convenient  than  M^g,  the  estimated  standard  deviation  is  2.5  times  as  large.  If 
the  Mg  index  for  two  rocks  differs  by  more  than  0.060,  there  is  a  68%  probability  that  the 
rocks  can  be  discriminated.  A  plot  of  Mg  versus  Rj  2  is  shown  in  Fig.  24.  The  Mg  parameter 
is  not  very  useful  for  acidic  rocks. 

5.2  GEOLOGIC  APPLICATIONS  FOR  THE  Rj  2  INFRARED  RATIO  IMAGING  METHOD 

The  two-channel  infrared  ratio  imaging  method  (described  in  previous  sections),  plus  a 
third  channel  for  temperature  correction,  is  a  new  tool  which  can  be  useful  for  discriminating 
among  silicate  rocks  and,  to  a  more  limited  degree,  for  identifying  unknown  silicate  targets 
as  belonging  to  a  certain  group  of  rocks.  Although  the  rock  groups  recognized  by  the  Rj  2 
ratio  are  different  from  groups  defined  by  traditional  rock  classification  schemes,  the  ratio 
does  considerably  reduce  the  number  of  rocktype  possibilities  which  an  unknown  target  can 
be.  Though  infrared  ratio  imaging  is  new,  there  are  some  geological  applications  for  which 
it  may  be  useful,  even  in  its  present  state  of  development. 

The  end  members  of  the  silicate  rocks,  namely  the  highly  acidic  and  ultrabasic  rocks, 
are  relatively  simple  to  discriminate  from  other  silicates.  Therefore,  the  most  straightfor- 
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TABLE  14.  REGRESSION  OF  THE  M,  MINERALOGICAL  INDEX 

D 

WITH  R  FOR  SET  C 


Mineral 

Parameter 

No. 

Coefficients 

Quartz 

1 

-.000798 

Anorthite 

2 

.000704 

Nepheline 

3 

.002673 

Olivine  (All) 

4 

.002264 

Magnetite 

5 

-.016910 

Orthoclase 

6 

-.000511 

Constant 

0 

1.208100 

Correlation  (r) 

.94 

Standard  Error  (s.e.) 

.019 
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ward  geological  applications  of  the  infrared  ratio  method  include  exploitation  of  these  end 
members.  At  one  extreme,  highly  acidic  rocks  can  be  identified  on  the  basis  of  their  low 
Rl,2  infrared  ratio,  such  as  the  quartz  sand  and  sandstone  shown  in  Fig.  25.  Three  obviously 
related  applications  are  exploration  for  construction  materials,  beach  surveys,  and  sand  dune 
mapping.  Sand  and  gravel  (chert  has  a  ratio  approximately  equal  to  that  of  quartz)  deposits 
can  be  easily  spotted  by  this  method,  in  night  or  day  operation.  Concerning  the  second  appli¬ 
cation.  the  difference  between  quartz  ard  calcite  ra'ios  is  so  pronounced  that  beach  surveys 
with  a  two-channel  thermal  infrared  scanner  should  permit  the  accurate  mapping  of  quartz- 
calcitc  sand  boundaries,  even  though  both  minerals  are  difficult  to  discriminate  by  eye  and 
practically  impossible  to  discriminate  by  photography.  Third,  quartz,  sand-dune  surveys 
likewise  should  he  aided  by  this  method  because  of  the  sharp  contrast  produced  in  R  be¬ 
tween  quartz  sand  and  most  other  natural  materials.  Quartz  sand  dunes  should  appear  distinct 
from  calcite  and  gypsum  sand  dunes,  also.  Al  the  other  extreme,  peridotites  and  other  ultra- 
basic  rocks  are  easily  discriminated  from  other  silicates  and  other  natural  materials,  a  fact 
which  leads  to  such  potential  uses  as  exploration  for  ophiolites  or  kimberlites.  Ophiolites, 
thought  to  he  remnants  of  fossil  subduction  zones,  are  basic  and  ul.rabasic  complexes  usually 
found  in  rugged,  mountainous  terrain,  where  remote  sensing  would  be  highly  desirable.  Kim¬ 
berlites,  which  may  be  economically  important,  should  also  he  identified  as  ultrabasic  com¬ 
plexes  by  the  infrared  ratio  method. 


There  are  types  of  applications  which  depend  primarily  on  the  infrared  ratio  for  discrimi¬ 
nation  alone  (not  for  identification).  One  is  the  mapping  of  volcanic  ash  flows  in  the  presence 
of  similarly  colored  alluvium.  Figure  26  is  a  blow-up  of  an  R]  g  rutio  image  from  a  region 
near  Pisgah  Crater,  California,  alongside  a  color  aerial  photo  taken  fimultaneously  with  the 
scanner  data.  The  rhyolitic  tuff  areas  (marked  T)  are  darker  than  the  less  acidic  alluvium 
in  the  infrared  ratio  image,  but  indistinguishable  from  the  alluvium  in  the  color  photo  Ash 
flows  have  traditionally  been  difficult  to  map  by  visible  means;  the  infrared  ratio  me.hixi  may 
be  helpful  in  solving  this  problem.  As  a  sidelight,  the  rhyolitic  tuff  in  this  intermittent  stream 
bed  near  Pisgah  Crater  contains  malachite  (a  copper  ore)  but  not  in  large  enough  quant, ties 
to  be  economically  important.  This  brings  up  the  subject  of  mineralogical  exploration,  which 
will  no  doubt  be  impacted  by  this  new  method  for  mose  problems  associated  with  great  local 
ffSi02  variations.  In  some  cases,  such  as  the  rhyolitic  tuff,  the  host  rock  for  the  ore  or 
precious  mineral  may  be  discriminable  from  other  rocks  in  the  scene  by  Rj  In  other  cases 
the  ore  itself  may  be  more  easily  discriminated  from  similarly  colored  rocks  by  the  infrared  ’ 
ratio.  For  example,  iron  formation  members  can  be  high  in  silica  content  as  compared  to 
dark  shales,  and  yet  the  two  can  be  equally  dark  in  the  visible  wavelength  region.  Concerning 
iron  and  several  other  metallic  ores,  the  R]  g  ratio  may  be  of  use  in  discriminating  widely 
different  n.etan, orphic  facies.  For  intermediate  temperatures  and  pressures,  the  green- 
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(a)  Channel  1:  8.2-10.9j.im 


lb)  Channel  2'  9.4-12.1 


FIGURE  25.  DISCRIMINATION  OF  ACIDIC  SILICATES 
NEAR  MILL  CREEK,  OKLAHOMA.  (Sand  quarry). 
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Infrared  Ratio  Image 


Aerial  Photo 


FIGURE  26.  COMPARISON  OF  THERMAL  INFRARED 
RATIO  IMAGE  WITH  AN  AERIAL  PHOTO  FOR  A  RHY¬ 
OLITIC  TUFF  (DARK  ON  RATIO  IMAGE,  MARKED  BY 
U  ASSOCIATED  WITH  MALACHITE.  Near  Pisgah  Cra¬ 
ter,  California.  North  is  toward  the  top  and  the  images 
are  approximately  1.5  km  on  a  side. 


i 
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set  ist  facies  contains  no  pyroxene  and  less  calcic  plagioclase  than  the  pyroxene  granulite 
facies,  which  does  contain  pyroxenes  [2G].  If  the  Mjg  mineralogical  index  (which  was  de¬ 
termined  for  igneous  silicates)  also  happens  to  lx1  valid  for  metamorphic  rocks,  the  green- 
schist  and  pyroxene  granulite  facies  may  lie  discriminate  by  the  ratio  imaging  method 

The  infrared  ratio  purely  as  a  discriminant  can  also  be  useful  in  general  geologic  map¬ 
ping,  especially  when  the  region  being  mapped  is  fairly  well  known  and  the  rock  types  are  of 
similar  color  —  but  quite  different  infrared  ratios  are  present  (e.g.,  some  anorthosites  and 
light-colored  felsite,  granite  and  limestone,  etc  ).  In  a  similar  way,  the  method  should  be 
useful  in  soil -mapping,  although  there  may  lie  arger  (SiC>2  -  Al^)  variations  with  depth 
of  soil  hori/.on  than  with  different  soil  species. 

Finally,  the  infrared  ratio  method  can  also  lx'  used  in  an  absolute  sense  to  help  identify 
silicate  rock  types  over  unknown  terrain.  For  instance,  if  a  rock  is  dark  in  the  visible  region 
and  has  an  Rj  2  ratio  near  1.2,  it  is  most  likely  a  basalt — and  definitely  not  a  peridotite. 
Using  the  visible  color  and  Rj  2  ratio  value,  the  rock-type  can  be  narrowed  to  a  few  possi¬ 
bilities.  This  would  be  most  useful  in  extraterrestrial  remote  sensing  and  in  reconnaissance 
mapping  of  inaccessible  terrain  on  earth  (where  ground  truth  is  sparse,  if  available  at  all). 
However,  since  atmospheric  effects  can  change  the  magnitude  of  the  infrared  ratio,  new  ra¬ 
tios  must  he  calculated  for  the  specific  experiment  performed.  In  remote  sensing  from  earth 
satellites,  for  instance,  the  total  earth  atmosphere  must  be  accounted  for:  this  problem  can 
and  will  be  attacked  both  by  extension  of  present  atmospheric  models  and  by  empirical  atmo¬ 
spheric  parameters  determined  from  Skvlab  infrared  spectrometer  measurements  [27  |. 

5.3  CONCLUSIONS 

A  method  has  been  presented  whereby  compositional  differences  among  silicate  rocks 
can  be  imaged  by  infrared  scanners  capable  of  recording  two  or  three  channels  in  the  8-  to 
14-pm  wavelength  region.  The  ratio  of  two  infrared  channels  has  been  related  to  chemical 
and  mineralogical  parameters  for  a  suite  of  widely  varied  igneous  silicate  rocks.  Infrared 
measurements  were  made  on  naturally  exposed  surfaces  of  these  rock  samples,  and  an  in¬ 
frared  ratio  (Rj  2)  was  calculated  for  each  rock  under  the  following  experimental  conditions: 
an  aircraft  altitude  1  km  above  ground,  a  clear  dry  atmosphere,  and  two  infrared  channels 
defined  by  a  particular  two-element  detector  currently  in  use. 


26.  W.  G.  Ernst,  Earth  Materials,  Prentice-Hall,  Inc.,  Englewood  Cliffs,  1969,  p.  135. 

27.  R.  K  Vincent,  The  NASA  Earth  Resources  Spectral  Information  System:  A  Data 
Compilation—  Second  Supplement,  Environmental  Research  Institute  of  Michigan,  Technical 
Report  No.  31650-156-T,  NASA  Contract  NAS9-9784,  1973. 
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Chemical  parameters,  consisting  of  oxide  weight  percentages,  were  measured  by  rapid 
rock  analysis.  Mineral jgical  parameters,  consisting  of  volume  percentages  of  natural  minerals, 
were  calculated  from  a  C.I.P.W.  computer  program,  with  the  oxides  as  inputs.  Forward  linear 
regressions  showed  that  chemical  parameters  are  better  correlated  with  Rj  than  mineralogical 
parameters  for  up  to  approximately  six  variables.  For  more  than  six  variables,  the  miner¬ 
alogical  parameters  are  better  correlated  with  Rj  r  The  best  single  parameter  is  %  SiO  ,  and 
the  best  two-parameter  combination  is  [SiC>2  +  5*4(Fe203)|,  with  (SiC>2  -  A120  )  a  close 
second  The  best  overall  compositional  parameter  is  a  mineralogical  index  (M  ),  which  is 
based  on  the  volume  percentages  of  1C  materials.  Another  mineralogical  index  (Mg),  which  is 
based  on  six  easily  identified  minerals,  was  created  for  use  with  rock  collections  which  have 
only  partial  modal  analysis  results. 

The  minerals  which  are  most  important  in  traditional  rock  classification  schemes  do  not 
have  the  same  relative  importance  In  relation  to  their  control  of  infrared  spectral  features. 

Thus,  the  R^  infrared  ratio  does  discriminate  among  silicate  rocks,  but  not  on  the  same 
basis  as  petrologists.  The  MJ6  mineralogical  index,  which  is  the  best  compositional  basis 
for  describing  the  Rj  2  discrimination,  was  divided  into  five  silicate  categories.  Only  three 
of  these  categories  are  comparable  to  traditional  rock  classifications. 

The  infrared  ratio  imaging  method  was  applied  to  two  test  areas  with  good  agreement 
between  scanner-measured  and  theoretically  calculated  values  of  R  In  the  area  around 
Mill  Creek,  Oklahoma,  quartz  sand  and  sandstone  were  readily  discriminated  from  other 
natural  materials  (including  other  silicates  such  as  topsoil,  and  nonsilicates  such  as  carbon¬ 
ates  and  vegetation).  In  the  area  around  Pisgah  Crater,  California,  dacites  and  basalts  were 
clearly  discriminated  and  were  Identifiable  on  the  basis  of  their  infrared  ratio  values.  On  a 
relative  basis,  a  rhyolitic  tuff  was  discriminable  from  surrounding  alluvium  on  the  basis  of 
the  R1>2  infrared  ratio,  whereas  the  two  could  not  be  separated  by  color  air  photographs  in 
the  visible  wavelength  region. 

Several  geologic  applications  have  been  suggested  for  the  infrared  ratio  method.  They 
include  the  following: 

(1)  exploration  for  construction  materials  (sand  and  gravel) 

(2)  surveying  beaches  (quartz  versus  calcite  sands)  and  mapping  sand  dunes 

(3)  exploration  for  ophiolites 

(4)  mineralogical  exploration 

(5)  mapping  of  volcanic  ash  flows 

(6)  reconnaissance  geologic  mapping  of  inaccessible  terrain 
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6 

SOIL  INFORMATION  FROM  REMOTE  SENSOR  IMAGERY 
6.1  INTRODUCTION 

Both  lame-scale  generalized  and  detailed  soil  information  are  required  for  tunnel  site 
selection  and  subsequent  tunnel  construction.  The  exact  location  and  method  of  detailed  soil 
investigations  are  determined  on  the  basis  of  preliminary  site  inspection  surveys  [28  |  and  in¬ 
terpretation  of  available  maps,  reports,  and  imagery.  It  is  the  purpose  of  this  section  to  de¬ 
scribe  some  of  the  types  of  soil  information  which  may  lx?  obtained  from  modern  remote  sensor 
data — particularly  multispectral  scanner  data. 

6.1.1  APPROACH 

Two  approaches  may  be  utilized  for  obtaining  soil  information  from  remote  sensor  data. 
One  approach  is  to  interpret  spatial  terrain  patterns  from  the  imagery  and  to  infer  the  nature 
and  extent  of  soils  on  the  basis  of  the  interpreter's  understanding  of  soil-environmental  re¬ 
lationships— particularly  those  of  vegetative,  topographic,  and  drainage  patterns.  This  is  the 
path  of  conventional  image  interpre  ation,  concerning  which  there  are  a  number  of  useful  ref¬ 
erences  [29-31 1 .  The  second  approach  is  to  analyze  only  those  portions  of  the  imagery  in 
which  the  soil  may  be  observed  directly.  Here  we  attempt  to  discern  soil  conditions  and  pat¬ 
terns  solely  from  the  appearance  of  these  bare  areas.  This  approach,  although  spatially  limited 
to  those  portions  of  the  terrain  in  which  soil  is  exposed  to  aerial  view,  is  particularly  suitable 
for  the  application  of  automated  processing  and  image  enhancement  techniques.  Consistent 
image  patterns  and  tones  are  assumed  to  represent  similar  soil  conditions;  these  may  then  be 
automatically  enhanced  or  selectively  delineated  with  current  computer  processing  techniques. 
But  since  these  techniques  are  designed  to  apply  to  specific  terrain  conditions,  their  effects 
on  non-bare  (vegetated)  areas  are  not  always  predictable.  Thus,  this  section  is  primarily 
concerned  with  lithological  applications  of  the  second  approach,  although  enhancement  of  vege¬ 
tative  differences  as  indicators  of  soil  conditions  is  also  discussed. 


28.  K.  Szechy,  The  Art  of  Tunneling,  Akademiai  Kiado,  Budapest,  1967. 


29. 


m  v  m’  FI?f\et  aL'’  Photointerpretation  of  Soils,  Manual  of  Photographic  Interpretation 
(R.  N.  Colwell,  ed.),  Am.  Society  of  Photogrammetry,  1960,  pp,  343-402. 


30.  V.  C.  Miller  and  S.  A. 
Surveys  and  Integrated  Studies' 
1968,  pp.  41-79. 


Schumm,  Aerial  Photographs  and  Surface  Features,  Aerial 
Proceedings  of  the  Toulouse  Conference,  UNESCO,  Paris, 


31.  A 
Films  for 
364. 


c”;  Clt7'b!rfma?11’  H  Gausman’  and  c-  L'  Wiegand,  Color  and  False  Color-IR 
Soils  Identification,  Photogrammetric  Engineering,  Vol.  37,  No.  4,  1971,  pp,  359. 
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6.1.2  BACKGROUND 

In  tunnel  site  selection,  soil  information  is  especially  important  for  shallow-depth  tunnels 
and  in  the  placement  and  construction  of  tunnel  portals.  Soil  information  may  also  be  impor¬ 
tant  as  an  indicator  of  subsurface  hydrologic  and  lithologic  conditions  over  the  lengtl  of  the 
tunnel  location.  In  other  words,  since  specific  soils  are  derived  from  certain  lithologic  parent 
materials  under  given  environmental  conditions,  natural  soil  permeability  and  drainage  may 
greatly  affect  subsurface  hydrologic  conditions.  In  general,  both  on-site  observations  and 
interpretation  and  analysis  of  aerial  images  are  required.  Precise  quantitative  methods  have 
been  developed  for  describing  and  sampling  soils  in  vertical  profiles  in  the  field.  Soil  samples 
are  also  extensively  analyzed  and  tested  in  the  laboratory.  The  challenge  is  to  extrapolate  these 
observations  to  the  surroundings  and  to  obtain  accurate  knowledge  of  important  spatial  varia¬ 
tions  in  soil  conditions  over  the  tunnel  site  location.  It  is  estimated  that  for  a  survey  of  inter¬ 
mediate  scale,  only  1  5,000,000  of  the  soil  is  actually  observed  from  borings  in  the  field.  Conse¬ 
quently,  the  soil  engineer  must  rely  on  the  art  of  image  interpretation  in  mapping  soils.  He 
not  only  requires  an  understanding  of  soil  morphology  ana  genesis,  but  more  importantly,  he 
also  must  know  the  relationships  of  soil  physics  to  geological  materials,  landforms,  topography, 
vegetation,  and  natural  drainage. 

Current  survey  practice  often  makes  use  of  aerial  photographs  as  a  mapping  base  and  for 
interpreting  the  area  over  which  field  observations  are  applicable.  In  some  initial  reconnais¬ 
sance  surveys,  aerial  photographs  largely  substitute  for  extensive  field  observations. 

In  working  with  images  we  have  recognized  certain  limitations  to  the  film-camera  system 
which  have  led  to  the  development  of  other  types  of  aerial  and  orbital  sensors.  Potentially  the 
newer  systems  provide  additional,  useful  terrain  imagery  to  the  soil  scientist  or  engineer. 

The  potential  of  one  such  system,  the  multispectral  scanner,  stems  from  the  large  amount  of 
spatial  and  spectral  data  collected  and  from  the  compatibility  of  these  data  with  modern  pro¬ 
cessing  facilities.  Some  of  the  information  concerning  soils  that  a  multispectral  scanner  sys¬ 
tem  may  prov'de  is  discussed  below.  In  addition,  a  practical  approach  for  realizing  this  po¬ 
tential  is  suggested. 

6.1.3  A  MULTISPECTRAL  SYSTEM 

For  purposes  of  this  discussion  it  is  important  to  briefly  consider  three  characteristics 
of  multispectral  systems:  the  spectral  range  and  resolution  over  which  they  are  sensitive,  the 
manner  in  which  the  data  are  recorded  and  stored,  and  several  of  the  methods  available  for  re¬ 
trieving  information  from  these  data. 
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Imagery  used  in  this  study  was  obtained  from  the  experimental  ERIM  M-7  airborne  scan¬ 
ner  system  (32  (.  This  system  synchronously  collects  12  bands  of  data  from  a  choice  of  19 
bands  over  the  0.32-  to  13.5-pm  spectral  range.  This  range  includes  radiation  in  the  ultra¬ 
violet,  visible,  near  infrared,  and  thermal  infrared  portions  of  the  electromagnetic  spectrum. 
The  scanner  has  a  maximum  2-milliradian  spatial  resolution  and  a  90°  total  field  of  view. 

The  12  bands  are  recorded  in  the  aircraft  on  compute  -compatible,  multitrack,  analog 
magnetic  tape;  this  allows  subsequent  electronic  processing  and  image  enhancement  of  these 
data.  Shown  in  Fig.  27  are  the  tape-recorded  signal  levels  associated  with  each  of  12  bands  for 
a  given  scanline  of  a  scanner  image.  Internal  radiation  reference  sources  which  allow  subse¬ 
quent  calibration  of  the  imagery  are  also  recorded  along  with  the  imagery.  The  data  illustrated 
were  recorded  from  the  Bureau  of  Mines,  multispectral-scanner  flight  over  a  portion  of  the 
Black  Hills  in  South  Dakota  (as  discussed  previously  in  this  report). 

The  simplest  method  of  retrieving  these  tape-recorded  data  is  to  photograph  sequentially 
a  CRT  video  image  of  each  of  the  12  spectral  Dands.  Each  image  portrays  the  terrain  in  a 
single  spectral  band  (as  shown  in  Fig.  27).  Contrasts  in  the  ultraviolet,  visible,  and  near- 
infrared  bands  are  a  result  of  solar  illumination,  terrain  aspect,  scanner  look -angle,  and  sur¬ 
face  reflectance  characteristics.  The  thermal  infrared  band  shows  contrasts  related  to  the 
absolute  temperature  and  emittance  characteristics  of  the  terrain  materials — warmer  areas 
are  characterized  by  lighter  tones.  Although  video  images  of  selected  bands  have  been  found 
useful  in  displaying  contrasts  associated  with  certain  terrain  features,  no  single  band  has 
proved  consistently  useful  for  identifying  soil  conditions.  However,  full  use  of  the  information 
contained  in  multispectral  data  requires  simultaneous  analysis  of  a  number  of  bands  in  combi¬ 
nation. 

Powerful  statistical  techniques  have  been  developed  for  use  with  computers  to  automatically 
recognize  scene  elements  on  the  basis  of  their  unique  spectral  characteristics  in  up  to  12  bands. 
These  techniques  have  been  used  experimentally  for  mapping  the  distribution  of  soils  in  areas 
bare  of  vegetation.  Several  limitations  of  this  technique,  however,  preclude  its  consideration 
as  an  operational  method  at  present.  These  include  the  large  amount  of  a  priori  soil  informa¬ 
tion  required  to  train  the  computer  and  validate  the  results  [33 1,  the  heterogeneity  of  terrain 


32.  P.  G.  Hasell,  Jr.,  Michigan  Experimental  Multispectral  Scanner  System,  Proceedings 
of  4th  Annual  Earth  Resources  Review,  NASA,  January  17-21,  1972,  pp.  34-1  through  34-13. 

33.  G.  Nagy,  G.  Shelton,  and  J.  Toloba,  Procedural  Questions  in  Signature  Analysis,  Pro¬ 
ceedings  of  the  7th  International  Symposium  on  Remote  Sensing  of  Environment,  Vol.  II,  May 
1971,  pp.  1387-1401. 
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(a)  IISGS  Topographic  Map  Section 


(b)  Ultraviolet  (0.33-0.38  pm) 


(c)  Violet  (0.41-0.48  pm) 


(d)  Blue  (0.4fi -0.49  pm) 

FIGURE  27.  MULTISPECTRAL  IMAGERY  DISPLAY.  Black  Hills 
South  Dakota,  21  May  1972,  1100  hr,  5000  ft  altitude. 
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(e)  Green  (0.50-0.54  jim) 
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FIGURE  27.  MULTISPECTRAL  IMAGERY  DISPLAY.  Black  Hills 
South  Dakota,  21  May  1972,  1100  hr,  5000  ft  altitude.  (Continued) 
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(i)  Infrared  (1.0-L4  jim) 


(j)  Infrared  (1.5 -1.8  4  m) 


(k)  Infrared  (2. 0-2. 6  Mm) 


(1)  Thermal  Infrared  (9.3-11.7  Mm) 


FIGURE  27.  MULTISPECTRAL  IMAGERY  DISPLAY.  Black  Hills 
South  Dakota,  21  May  1972,  1100  hr,  5000  ft  altitude.  (Concluded) 
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surface  which  limits  the  extension  of  statistical  recognition  criteria  [34  |,  and  the  large  costs 
associated  with  data  analysis  and  recognition  processing.  These  are  discussed  further  in 
Section  f  .3 

For  tunnel-site  selection  purposes,  we  require  a  procedure  that  produces  spatial  soil  in¬ 
formation  which  is  not  readily  available  from  less  expensive  aerial  photographic  systems, 
allows  use  of  the  engineer's  knowledge  of  soils  and  local  terrain  conditions  prior  to  and  subse¬ 
quent  to  image  production,  and  is  economically  cost  effective.  Thus,  Section  6.4  of  this  report 
is  concerned  with  the  constraints  on  operational  multispectral  systems  for  identifying  and 
mapping  soil  conditions. 

6.1.4  SOIL  REFLECTANCE 

Implicit  in  the  use  of  imaging  systems  for  identification  of  soil  conditions  is  the  assumption 
that  the  reflectance  or  emittance  characteristics  of  soils  are  wavelength-dependent  and  rela¬ 
tively  specific  for  differing  physical  conditions.  In  other  words,  the  physical  and  chemical 
nature  of  soils  affect  certain  radiation  absorption  characteristics,  and  these  effects  are  ob¬ 
servable  from  aerial  platforms. 

In  selecting  spectral  bands  to  use,  we  must  know  something  about  the  spectral  character¬ 
istics  of  soils.  Very  general  relationships  of  soil  conditions  to  their  reflectances  have  long 
been  noted  but  seldom  quantitatively  investigated.  Soil  reflectance  generally  increases  mono- 
tonically  with  increasing  wavelength  over  the  ultraviolet,  visible,  and  near  infrared  portions  of 
the  spectrum.  Most  frequent  observations  concern  the  more-or-less  uniform  decrease  in  soil 
reflectance  with  increasing  moisture  content,  increasing  percentages  of  organic  matter,  and 
increasing  particle  or  aggregate  sizes.  Also  noted  are  differences  in  visual  color  related  to 
the  form  and  amount  of  iron  oxides. 

In  a  study  of  160  surface  soil  samples  collected  from  36  states,  Condit  [35 1  showed  only 
three  significant  variations  in  the  shapes  of  spectral  reflectances  curves  for  the  0.3-  to  1.0- 
/an  range:  a  relatively  linear  increase  with  increasing  wavelength  (Chernozen  type),  a  convex 
shape  resulting  from  lower  rates  of  increase  in  reflectance  beyond  about  0.6  (im  (pedalfer 
type),  and  an  inverted  S-shaped  curve  resulting  from  a  leveling  or  slight  decrease  and  then  an 
increase  in  reflectance  beyond  0.7  /an  (lateritic  type).  In  other  words,  only  slight  variations 


34.  R.  F.  Nalepka  and  J.  P.  Morgenstern,  Signature  Extension:  An  Approach  to  Opera¬ 
tional  Multispectral  Surveys,  Environmental  Research  Institute  of  Michigan,  Report  No.  3165- 
152-T,  1973. 

35.  H.  R.  Condit,  Application  of  Characteristics  Vector  Analysis  to  the  Spectral  Energy 
Distribution  of  Daylight  and  Spectral  Reflectance  of  American  Soils,  Applied  Optics,  Vol.  11, 
No.  1,  1972,  pp.  74-86. 
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in  the  shapes  of  reflectance  curves  appear  to  occur  as  a  result  of  differences  in  soil  conditions 
Average  reflectance  (albedo)  levels  vary  quite  greatly,  however. 

Color  contrasts  have  long  been  used  as  an  easily  observable  soil  property  which  may  indi¬ 
cate  more  important  soil  differences.  Differences  in  surface  color  may  indicate  the  locction 
of  different  soil  types  or  degrees  of  surface  erosion.  (Some  of  the  older  soil  taxonomic  sys¬ 
tems  maintain  a  strong  emphasis  on  color  as  a  classification  criterion.)  Since  the  late  1940's, 
soil  colors  have  been  routinely  recorded  in  the  field  through  use  of  the  serai-quantitative 
Munsell  color  standards  [36  ] . 

6.2  RATIOED  IMAGES 

Intermediate  to  the  simple  playback  of  single  video  images  and  the  multiband  statistical 
recognition  of  scene  features  is  the  ratioing  of  two  bands  of  multispectral  data.  Using  an  image 
processor,  it  is  a  relatively  inexpensive  and  simple  procedure  that  divides  ,h<.  recorded  sig¬ 
nals  in  one  band  by  those  in  another  to  produce  a  new  (ratioed)  image. 

This  new  image,  which  is  a  derivative  of  the  two  original  bands,  displays,  by  variations 
in  tone,  those  areas  which  have  differing  amounts  of  spectral  change  from  one  band  to  the 
other.  In  other  words,  the  contrasts  inherent  in  the  ratioed  image  result  from  the  relative 
differences  between  the  two  bands  and  are  not  dependent  on  absolute  signal  levels  in  either 
band.  This  fact  is  important.  A  ratio  image  is  fundamentally  different  from  a  single-channel 
video  image,  in  spite  of  the  superficial  similarity  of  the  ratioed  image  to  a  single-band  video 
image — or  to  an  aerial  photograph  for  that  matter.  The  ratioed  image  is  less  dependent  on 
illumination,  temperature,  and  directional  variations  of  the  scene,  and  it  is  more  nearly  a 
function  of  terrain  reflectance  or  emittance  characteristics  than  the  radiance  levels  recorded 
in  any  one  band.  The  ratio  also  tends  to  emphasize  subtle  spectral  differences  when  signal 
levels  in  the  denominator  are  low — a  feature  of  the  hyperbolic  nature  of  ratios. 

6.2.1  MUNSELL  HUE 

In  an  experiment  (for  the  USDA  Experiment  Station  at  Weslaco,  Texas)  to  delineate  soils 
from  multispectral  data,  it  was  noted  that  the  11  soil  types  differed  in  their  redness,  or  domi¬ 
nant  Munsell  hue  [37 1.  Eight  were  in  the  5  YR  (Yellow-Red)  range  and  three  were  in  the  10  YR 
range. 


36.  J.  A.  Shields,  et  al.,  Measurements  of  Soil  Color,  Canadian  Journal  of  Soil  Science, 
Vol.  46,  1966,  pp.  83-90. 

37.  A.  H.  Munsell,  A  Color  Notation,  Munsell  Color  Co.,  Inc.,  Baltimore,  1947. 
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Two  types  of  spectral  curves  were  obtained  from  multispectral  data  of  20  soil  sample  lo¬ 
cations.  Illustrated  in  Fig  28,  these  two  types  of  curves  clearly  distinguish  the  10  YR  from 
the  l  edder  i  YU  soils,  although  considerable  variation  in  the  absolute  levels  occurs  within  each 
t\pc.  (The  maximum  and  minimum  for  each  tyjx  are  shown.)  The  primary  difference  between 
the  two  types  of  curves  occurs  in  the  rate  of  increase  between  the  0.52  pm  (green)  and  the  0.62 
I  m  (red)  wavelength  range.  The  redder  5  YR  soils  more  than  double  their  reflectance,  while 
the  10  YR  soils  remain  about  the  same  over  this  range.  This  fact  suggested  that  the  soils  of 
this  area  could  lx  separated  on  the  basis  of  Munsell  hue  by  a  simple  ratio  of  the  0.61  to  0.70 
/  m  (red)  band  to  the  0,50  to  0.54  pm  (green)  band.  Figure  29  compares  an  annotated  video 
image  showing  the  locations  of  the  soils  with  the  ratioed  image  of  these  two  bands.  The  darker 
soil  areas  correspond  to  those  with  10  YR  hues,  while  the  lighter  areas  have  5  YR  hues  on  the 
ratioed  image.  It  is  likely  that  the  differences  in  hue  result  from  the  relative  ferric  iron  con¬ 
tent  of  the  soil  parent  materials,  but  the  sample  analysis  is  not  yet  complete  enough  to  confirm 
this  hypothesis.  Obvious  differences  related  to  the  average  lightness  or  darkness  of  these  soils 
(Munsell  value)  were  reduced  with  this  ratio  image.  Contrasts  related  to  Munsell  value, which 
are  apparent  on  the  single-band  image,  are  probaMy  related  to  organic  matter  and  relative 
moisture  contents. 

6.2.2  NATURAL  SOIL  DRAINAGE 

Poorly  drained  soil  areas  generally  appear  darker  than  surrounding  drier  soils  in  the 
reflectance  wavelengths  /Jin).  During  the  day,  this  contrast  is  often  reversed  in  the  thermal 
infrared  band  that  is,  moist  areas  containing  more  organic  matter  appear  warmer  (because 
of  their  darker  color)  than  drier  soil  areas.  A  ratio  of  a  reflective  band  and  a  thermal  band 
then  serves  to  enhance  differences  associated  with  soil  moisture  and  organic  matter  accumula¬ 
tion  usually  an  indication  of  natural  soil  drainage.  Shown  in  Fig.  30  is  the  coincidence  of  a 
dark  pattern  in  a  bare  field  with  a  soil  mapped  as  poorly  drained.  Differences  between  well 
drained  and  somewhat  poorly  drained  are  not  so  obvious,  however. 

6.2.3  VEGETATION  VIGOR 

The  ability  of  a  ratioed  image  to  enhance  tone  differences  associated  with  terrain  conditions 
depends  on  the  existence  of  relative  spectral  differences  in  two  or  more  bands.  The  most  effec¬ 
tive  differences  are  tone  reversals  from  one  band  to  another.  For  example,  comparisons  of 
vegetation  in  the  0.62-  to  0.70-pm  (red)  chlorophyll -absorption  band  and  the  0.70-  to  0  94 -pm 
(near-infrared)  reflectance  band  show  that  vigorous  vegetation  is  darker  in  the  red  band  and 
lighter  in  the  infrared  band  than  stressed  or  less-dense  vegetation.  Examples  of  these  dif¬ 
ferences  in  lawn  vigor  are  shown  for  a  small  tcwn  in  western  Virginia  in  April,  1971.  Figure 
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Thermal  Image 
9.3-11.7  ,i  m 


Near-IR  Image 
1.5-1. 8  pm 


Ratio  Image 
9.3-11.7  pm 
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FIGURE  30.  ENHANCEMENT  OF  NATURAL  SOIL  DRAINAGE 
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31  compares  l,laek-i„„l-wl,i,o  aerial  pi, . .  will,  a  noa,-m  rod  raliocd  imago  ,llcsc 

,u'r‘'  |,,'"‘"*«>'d  '»'■  l"lc:,l  Highway  Administration  |3«|. 

fi.2.^  SCENE  CONTRAST  SUPPRESSION 

Aerial  I, paces  l.vp.cally  slum  terrain  eopirasts  related  larcoly  lo  land  use  practices  Ip 

acricullural  areas.  hare  Me], Is  otlea  are  . .  cop.parisoa  fields  win,  dark  cree„ 

ceops  These  sharp  eoplrasls  have  II, e  etfee,  ot  conceal!, ip.  olher.  more  sahue,  lerrain  pal¬ 
let,,,  related  lo  drai„a«e.  physiography,  or  „a,„ral  vece,alio„  dislribulion.  Also,  II, e  occur- 
mu  i  of  intensive  dark  areas  on  the  imagery  makes  it  less  suitable  as  a  field  mapping  base 
hecause  of  II, e  ilifflcplly  ip  sceinc  l*n  or  imnil  a.molallon.  In  those  dark  areas. 

Flit uro  32  demons, rale,  thnt  II, e  ahilily  oxisls  reduce  these  marked  contrasts  sicnificanllv 
iniaces  of  two  spectrally  similar  liands  (green  and  vjolel  in  IIP,  case,  Usingdif- 
feren,  hands  has  the  effect  of  emphasisinc  or  suppress!,,,,  different  comhinalions  o,  terrain 
patte,  ns  different  degrees.  The  ralioed  figure  In  this  example  reduces  most  „,e  contrast 
Wilhm  the  cultivated  areas  hut  luainlains  contras,  for  Ihe  wooded  porlions  o,  Ihc  Image  The 
'.die,  areas  are  scarcely  discernible  on  Ihe  single-channel  video  image-a  panchroma, ic 
phot ogiapli  would  appear  similar  to  this  video  image. 

6.2.5  SOIL  EMITTANCE 

Energy  radiated  from  a  soil  surface  in  the  thermal  infrared  portion  of  the  spectrum  (^3 

,  'S  Pr,manIy  8  fUnC"0n  °f  ",0  emlltanc<‘  *'>*»'**  temperature  of  the  soil.  Absolute 
ompera'ure,  in  turn,  is  dependent  on  the  thermal  properties  of  the  soil  and  the  diurnallv 
changing,  local  climatic  conditions. 

Lniillance  charnolcrislics  of  soils  (II, a,  |s,  how  closely  they  simnlale  a  Iheorelically  per- 
eu  emiller  a  Waekhody")  are  wavelengll, -dependent.  In  Ihe  »-  lo  10-,,,,,  range,  soil  Ihormal 
emillnnce  ,s  lowered  l,y  II, e  so-called  reslslrahlen  eftecls  o,  silicale  minerals.  TIPs  lower 
emiltancc  reduces  Ihe  apparen,  radian,  energy  (hlackhody  lemperalu,  o)  in  this  range  In  com- 
pari  son  lo  II, e  apparen,  radian,  energy  a,  on, or  wavelengths  for  „,0  same  soil.  Field  observa¬ 
tions  rcporled  a,  Purdue  |39|  show  Ihe  decrease  in  apparent  lempernlnre  in  Iho  «-  10-,.,,', 

Muge  fo.  mineral  sous  and  a  sandy-gravel  road  (see  Fig.  331.  Grealer  apparen!  temperature 


Malerhals^a/Virgini'a  Tesl'siles  ' r'  ".'""’e  W<'nUncali™  °>  Terrain  Features  and 
search  hislilule  of  Michigan,  Itcporl  No.  27600-9™  ITO°'  Tecl,0,(|Ul'“'  Environmental  Re- 

197 1 3^  Al'nU'“'  N°-  "•  Talioralory  for  Agricullural  Remole  Sensing,  Purdue  Umycrsily, 
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\erjm 


0.41-0.48  /jm 
Video 


0.50-0.54  (Green) 
0.41-0.48  Violet) 


FIGURE  32.  COMPARISON  OF  A  VIDEO  IMAGE  OF  AGFilCULTURAL 
LAND  USE  WITH  A  RATIOED  IMAGE  OF  THE  SAME  AREA.  Ratio 
image  is  designed  to  suppress  scene  contrast  associated  with  cropped 
versus  bare  fields.  Dark  area  in  upper  right  is  a  cloud  shadow. 


101 


ORMERl 


WILLOW  RUN  LAftOWATOkllS  THE  UNIVERSITY  Of  MIC  MIGA 


WAVELENGTH 


FIGURE  33.  BLACKBODY  TEMPERATURES 
SHOWING  THE  RESTRAHLEN  EFFECT  FOR 
MINERAL  SOILS  IN  THE  8  TO  10  4 m  WAVE¬ 
LENGTH  RANGE.  Data  from  LARS  (Purdue)  — 
annual  report  No.  4. 
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differences  are  indicated  for  the  two  silt  loam  soils  than  for  the  clay  loan,  soil,  and  the  organic 
muck  s(Ml  sh(nvs  litt,e  variation  in  the  reststrahlen  band.  Also,  the  crop  vegetation,  hard- 
surfaced  road,  and  water  surface  show  little  or  no  apparent  temperature  variation  in  this 
range  Absolute  temperature  values  are  not  of  concern  here. 

Roughly,  the  greater  the  SiO.,  content  of  a  soil  surface,  the  lower  the  soil  omittance  is 
expected  to  Ik  ,n  the  H-  to  10-/.  m  range.  Since  SiO.,  frequently  makes  up  the  bulk  of  the  sand- 
s,/c  fraction  of  well -weathered  soils,  a  rough  correlation  between  a  ratio  of  two  thermal  bands 
and  surface  soil  texture  may  be  possible.  Indeed,  the  reststrahlen  effect  has  been  shown  to 
vary  somewhat  with  particle  size  as  well  as  with  mineral  content  [40], 

To  test  this  approach  to  remote  determination  of  soil  texture,  two  bands  of  thermal  infra¬ 
red  imagery  collected  near  Weslaco,  Texas,  were  ratiood.  The  two  bands  used  wore  in  the 
H  '  to  10.9-/  m  and  the  9.4-  to  12.1-pn.  wavelength  ranges  -the  former  being  the  reststrahlen 
hand  for  silicate  minerals.  Results  are  presented  in  Fig.  34,  Four  soil  types  ranging  from 
clav  loam  to  loamy  fine  sand,  plus  field  crops,  are  shown  in  the  scene.  W'hile  it  would  not  be 
possible  to  identify  these  soil  types  individually  from  the  single  thermal  images,  it  appears 
that  the  ratiood  image  does  indicate  by  tone  contrast  (darker  areas  have  sandier  textures) 
the  locations  of  differing  soil  types.  Soil  samples  collected  from  this  area  are  currently  being 
analyzed  and  should  provide  confirmation  of  this  conclusion.  It  should  bo  noted  that  differences 
in  soil  moisture  and  in  cultivation  practices  are  likely  to  affect  radioed  image  results  — in  addi¬ 
tion  to  any  c‘ff<-,ct  produced  hy  the  silicate  composition  of  these  soils. 

G.3  MULTICHANNEL  RECOGNITION 

If  soils  or  other  materials  of  interest  are  sufficiently  different  from  one  another  in  their 
spectra,  they  should  be  discnniinablc  in  some  manner.  Consider  two  objects,  A  and  B,  in  two 
spectral  bands  (dimensions),  as  shown  in  Fig.  35.  These  objects  of  interest  may  be  recognized 
by  ascertaining  that  a  sample  occurs  within  some  boundary  about  the  objects  in  the  space  (/, 

'V-  1,1  tl,is  sample,  little  confusion  is  possible  since  the  objects  do  not  overlap.  They  may 
be  discriminated  by  very  simple  boundaries,  such  as  a  ratio  value.  But  more  complex  situa¬ 
tions  arise  in  nature:  here,  simple  two-dimensional  boundaries  can  cause  misclassification  of 
scene  objects  of  interest.  For  example,  with  ratio  discrimination,  note  that  portions  of  object 
C  would  be  mistakenly  classified  as  either  object  A  or  B.  In  such  cases,  a  more  general  type  of 
boundary  is  necessary,  usually  based  on  a  statistical  model  in  which  the  probabilities  orlikeli- 


nf  ,  40'  R  i  Ly°n’  Paluation  of  I,lfrar0<i  Spectrophotometry  for  Compositional  Analysis 
of  Lunar  and  Planetary  Soils,  NASA  Report  No.  49104,  Stanford  Research  Institute,  1964 
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8.2  to  10.9  it  m 
and  9.4  to  12.1  ft m 
(ratio) 


Area: 

A  -  Rio  flay  loam 
B  -  Ramadero  sandy  clay  loam 
C  -  Hidalgo  fine  sandy  loam 
D  -  Hebbronville  loamy  fine  sand  (?) 
E  -  Crop  vegetation 


FIGl'RE  34.  COMPARISON  OF  TWO  THERMAL  IMAGES  WITH 
A  RATIOEr  IMAGE  OF  A  SOIL  TEST  AREA  NEAR  WESLACO. 

TEX  vS 
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floods  of  the  various  objects  are  compared  for  each  sample  to  be  classified.  It  is  often  neces¬ 
sary  that  tin  spectral  characteristics  nt  several  objects  be  considered  simultaneously  to  obtain 
the  best  results.  Multi-channel  analog  and  digital  computer  techniques  have  been  developed  for 
tin  s  purpose 

In  situations  where  the  spectral  distributions  of  objects  overlap  or  are  expected  to  overlap 
(generally  the  case),  simple  ratio  techniques  do  rot  pi  ovtde  adequate  discrimination  capabilities. 
The  concept  of  likelihood  making  a  decision  on  the  basis  of  probability  —  is  more  appropriate. 
If  the  two  distributions,  A  and  B,  overlap  (see  Vig.  36),  equi-probable  contours  are  drawn  for 
the  distributions  of  objects  A  and  B  At  the  intersections  of  those  contours  of  equal  probability 
value,  A  and  B  are  equally  likely  or  probable:  and  a  curve  may  be  drawn  through  these  points 
to  form  a  decision  surface  on  either  side  of  which  decisions  mav  be  made  for  any  point.  The 
function  used  for  this  decision  is  known  as  the  likelihood  ratio: 


L(A  A  j,  a2) 


|)(A  A  j ,  a2) 
p(B  A  j,  A 2 ) 


(ID 


where  L  the  likelihood  ratio 

pfA  A  j,  a2>  die  probability  of  object  A,  given  measurements  (Aj,  Ag) 
p(B  Aj,  a2)  the  probability  of  object  B,  given  measurements  (A  , 

If  the  decision  rule  is  established  sucli  that 


p(A  Aj,  A2) 

1j  “  p(BlAj,  A2)  '  1 


(12) 


tlien  this  is  a  likelihood-ratio  decision,  and  the  bounds  are  represented  by  this  decision  rule. 
Fast  computers  allow  the  extension  of  this  likelihood  decision  concept  to  a  number  of  dimen¬ 
sions  (spectral  bands)  and  multiples  of  objects. 

Using  soil  sample  areas  selected  from  southeastern  Michigan,  an  analog  image  processor, 
SPARC,  was  used  to  implement  this  likelihood-ralio  decision  rule.  The  processor  accepted 
data  from  four  spectral  bands  for  three  surface  soils  and  applied  the  likelihood  test  to  deter¬ 
mine  the  distribution  of  these  materials  for  the  flightline.  The  assumption  is  made  that  the 
probability  distribution  of  the  samples  is  approximated  by  a  normal  or  Gaussian  distribution. 

Signals  from  each  of  the  four  selected  spectral  bands  are  supplied  to  summing  amplifiers 
from  a  tape-loop  recorder,  along  with  the  mean  values  for  the  particular  object  of  interest. 
The  level  at  the  output  of  the  first  summer  is  set  to  a  particular  value  by  manually  adjusting 
a  potentiometer  (1  jdj).  Outputs  from  the  summers  are  repetitively  examined  on  an  x-y  os¬ 
cilloscope  an  i  adjusted  to  show  no  correlation.  The  level  of  the  signal  at  the  output  of  the 
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second  summer  is  set  to  the  same  value  as  the  output  of  the  first  by  adjusting  the  second  po¬ 
tentiometer,  and  so  on.  These  signals,  now  uneorrelated,  are  squared  and  added  to  obtain  the 
exponent  of  the  normal  distribution.  The  output  is  then  exponentiated  to  obtain  the  probability 
of  all  the  signals. 

The  same  operations  are  performed  on  the  signals  for  all  other  sample  material  distribu¬ 
tions.  The  outputs  are  then  supplied  to  the  final  summer  with  the  proper  sign  and  tested  with 
a  threshold  device.  If  the  unknown  output  signal  is  greater  than  a  present  threshold,  an  object 
of  interest  is  registered  as  present.  This  final  output  signal  is  displayed  on  a  CRT  and  pho¬ 
tographed  with  a  filmstrip  printer.  Thus,  the  output  is  a  black  (object)  and  clear  (not  object) 
recognition  filmstrip  showing,  point  by  point,  the  results  of  the  likelihood  test  for  each  ma¬ 
terial  in  the  entire  study  area. 

Two  examples  of  multichannel  recognition  of  soil  areas  are  included  here  for  illustration 
both  areas  are  in  southern  Michigan.  Figure  37  shows  recognition  of  three  soil  types  having 
different  drainage  characteristics.  For  purposes  of  comparison,  a  drainage  map  (based  on  ex¬ 
isting  soil  survey  maps)  is  provided.  The  computer-generated  distribution  of  well-drained 
and  somewhat  poorly  drained  soils  compares  well  with  the  soil  drainage  map,  although  there 
is  considerably  greater  detail  in  the  processed  imagery.  The  somewhat  poorly  drained  soils 
occur  on  the  lower  slopes  of  this  undulating  terrain.  Poorly  drained  soils  recognized  by  the 
computer  constitute  too  small  an  area  to  >e  differentiated  on  the  conventional  map.  Field  ex¬ 
amination  determined  that  these  areas  were  correctly  delineated.  The  multispectral  data  used 
to  produce  these  images  were  collected  at  2000  ft  altitude.  The  following  spectral  bands  were 
employed:  0.48-0.52  pm  (blue-green),  0.67-0.94  pm  (near-IR),  2. 0-2. 6  pm  (mid-IR),  and  9.3- 
11.7  /mi  (thermal-IR). 

The  currently  orbiting  Earth  Resources  Technology  Satellite  (ERTS)  also  collects  multi- 
spectral  data,  but  at  a  much  coarser  resolution  because  of  its  500  mi  altitude.  Figure  38  shows 
an  attempt  to  differentiate  different  soil  drainage  classes  by  using  ERTS  multispectral  data. 

In  this  example,  correlation  of  the  digital  computer-generated  results  with  existing  soil  maps 
indicates  that  consistently  correct  identification  of  these  bare  fields  occurred  only  in  the  case 
of  the  mineral  versus  organic  (peat  and  muck)  soils.  These  data  were  collected  at  a  sub¬ 
optimum  time  of  year  for  soils  differentiation  (less  than  10^  of  the  ground  was  free  of  vege¬ 
tative  cover,  as  indicated  by  sparse  soil  recognition).  Research  is  continuing  into  the  use  of 
orbital  multispectral  imagery  to  separate  soil  drainage  classes. 

6.4  TOWARD  AN  OPERATIONAL  SYSTEM 

Earlier,  we  suggested  several  requirements  which  should  be  met  if  we  are  to  consider 
multispectral  systems  for  an  operational  role  in  soil  surveys.  Ultimately,  we  must  be  con- 
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Soil  Drainage  Map  Well  Drained 

(Miami  Series) 


Somewhat  Poorly  Drained 
(Conover  Series) 


FIGURE  37.  COMPARISON  OF  A  CONVENTIONAL  SOIL  DRAINAGE  MAP 
WITH  SOIL  RECOGNITION  IMAGES  FOR  FLEMING  CREEK  WATERSHED 
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it  mod  with  such  applications  if  wo  arc  to  justify  continued  public  support  for  the  development 
of  sucli  systems.  What  are  some  of  the  limitations  in  the  applic:  m  o.  nuiltispectral  data  for 
identifying  and  mapping  soil  conditions'’ 

One  limitation  is  that  often  only  a  portion  of  the  soil  is  exposed  to  aerial  view  in  an  area 
to  be  surveyed  To  maximize  bare  soil  areas,  it  is  necessary  to  collect  remote  sensor  data 
from  agricultural  areas  either  in  the  fall  or,  preferably,  in  the  spring  prior  to  substantial  crop 
growth  Even  then,  there  will  he  areas  of  perennial  vegetation  or  cultural  development  for 
which  no  direct  soil  information  can  be  obtained  from  aerial  imagery.  Aerial  imagery  can 
delineate  only  surface  or  near  surface  conditions:  taxonomic  soil  units  established  from  sub¬ 
surface  profiles  cannot  be  identified  from  imagery.  Precise  identification  of  soils  must  still  be 
made  on  the  ground.  Remote  sensors  can  help  extrapolate  those  observations,  reduce  their 
numbers,  and  provide  easier  determinations  of  soil  boundaries;  but  they  cannot  replace  the 
task  of  field  sampling. 

The  processed  results  presented  here  were  produced  in  a  research  mode  using  an  experi¬ 
mental  scanner  and  playback  system.  Other  techniques  may  prove  to  lie  better  suited  for  en¬ 
hancement  of  these  or  other,  more  useful  soil  conditions.  For  example,  a  ratio  of  two  fre¬ 
quencies  of  radar  data  hopefully  may  provide  new  synoptic  information  concerning  near-surlam 
properties  of  different  soils.  Experiments  are  currently  underway  to  test  this  idea. 

hi  time,  we  hope  to  attach  numerical  ranges  to  absolute  values  for  certain  soil  conditions. 

If  successful,  this  may  make  possible  an  automatic  discrimination  process  to  select  and  print 
out  recognition  images  of  soil  conditions  from  selected  ratios.  Here,  discrimination  could  be 
accomplished  automatically  by  a  simple  level -slicing,  or  thresholding,  procedure. 

6.4.1  SPECTRAL  VARIATION’S 

In  addition  to  enhancing  contrasts  associated  with  soil  conditions,  ratioed  images  reduce 
certain  scene  variations  which  occur  as  a  result  of  illumination  or  directional  reflectance 
differences.  Recent  studies  have  shown  that  soils  are  not  uniformly  diffusing  Lambertian 
surfaces.  When  the  sun  angle  is  low,  the  same  soil  can  have  quite  different  appearances 
depending  upon  the  direction  of  view  (see  Fig.  39).  Variations  in  terrain  slope  and  aspect  can 
have  a  similar  effect  even  when  the  sun  is  high.  These  effects,  associated  with  solar  direction, 
look  angle  geometry,  and  terrain  slope,  are  recorded  by  the  scanner.  Frequently,  one  side  of 
a  video  image  is  uniformly  lighter  than  the  other  side  as  a  result  of  scanning  alternately  into 
and  away  from  the  sun’s  direction.  These  scan  angle  variations  are  only  slightly  wavelength- 
dependent  .  A  ratio  of  two  bands,  particularly  two  bands  which  are  spectrally  close,  effectively 
i  luces  this  directional  reflectance  variation,  Ratios  of  adjacent  bands  have  been  used  for  a 

number  of  years  as  a  preprocessing  technique  to  eliminate  gross  directional-reflectance  varia¬ 
tions. 
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Solar  illumination  may  also  vary  as  the  aircraft  collects  data  along  a  flightline.  Again, 
these  variations  are  not  strongly  wavelength -dependent,  and  a  ratio  of  two  reflective  bands  helps 
eliminate  cloud  shadows  and  other  types  of  illumination  variations. 

6  4  2  CONCLUSIONS 

We  have  attempted  to  present  some  of  the  kinds  of  soil  information  available  from  selec¬ 
tively  enhanced  multispectral  images  —  soil  color,  terrain  patterns,  natural  drainage,  and  soil 
texture.  The  relative  simplicity  of  the  ratio  approach  seems  to  commend  it  for  study  as  a 
source  of  soil  survey  imagery.  On  the  other  hand,  experience  with  single  video  images  indicates 
no  more  soil  information  than  is  contained  in  aerial  photography.  Statistical  recognition  pro¬ 
cessing  is  laborious,  and  additional  research  into  its  use  in  standard  survey  practice  is  re¬ 
quired. 

Itatioed  images  are  produced  at  the  same  real-time  rate  as  single-band  video  playbacks 
and  at  a  comparable  cost.  If  the  soil  surveyor  can  suggest  what  kinds  of  surface  soil  and  ter¬ 
rain  conditions  would  he  useful  for  discriminating  and  mapping  soils  in  his  area,  he  might  be 
able  to  specify  the  ratios  most  likely  to  emphasize  any  soil  differences.  When  the  resulting 
enhanced  images  are  |  reduced,  he  can  then  use  the  full  range  of  his  knowledge  and  experience 
in  interpreting  them.  Other  ratios  of  multispectral  bands  used  for  enhancing  vegetation  dif¬ 
ferences  have  been  reported  recently  [41,  42  | .  These  may  also  prove  useful  to  the  soil  sci¬ 
entist  in  terrain  areas  where  he  cannot  observe  the  soil  surface  directly. 


41  R.  L.  Pearson  and  L.  D.  Miller,  Remote  Mapping  of  Standing  Crop  Biomass  for  Esti¬ 
mation  of  the  Productivity  of  the  Shortgrass  Prairie,  Pawnee  National  Grasslands,  Colorado, 
Proceedings  of  the  8th  International  Symposium  on  Remote  Sensing  of  Environment  Vol.  II, 
October  1972,  pp.  1355-1379. 

42.  C.  I  Olson,  Jr.,  Remote  Sensing  of  Fomes  annosus  in  Forest  Stands,  Proceedings  of 
the  8th  Inter  lational  Symposium  on  Remote  Sensing  of  Environment,  Vol.  II,  October  1972,  pp. 
1381-1384. 
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GENERAL  CONC LESIONS  AND  RECOMMENDATIONS 


Remote1  sousing  for  geologic  recomiaissiince  is  on  t he  verge  of  Incoming  operational.  ERTS-1 
lias  already  pt  r  ormed  the  two  tasks  of  alerting  tin1  user  community  to  remote  sensing  and  pro¬ 
viding  data  with  winch  to  develop  useful  analysis  procedures.  Because  of  the  wide  publicity  given 
the  successful  experiments  of  ERTS-I,  the  authors  believe  that  any  further  justification  for,  de¬ 


scription  of,  or  proof  of  the  usefulness  of  remote  sensing  for  geologic  purposes  would  be  redun¬ 
dant.  Therefore,  we  have  concentrated  on  contributing  directly  to  the  further  development  of 
geologic  remote  sensing.  To  this  end,  a  test  site  was  selected  in  the  Black  Hills  of  South  Dakota, 
and  MSS  airborne  data  were  collected  during  two  flights  in  May,  1972.  Although  the  cost  of  gatli- 
ei  mg  complementary  SLAB  data  expressly  tor  this  contract  was  prohibitive,  Strategic  Air 
Command  SEAR  imagery  was  obtained  for  a  portion  of  the  area  flown:  it  could  not  be  included 
in  this  report  because  of  analysis-time  limitations. 


We  feel  that  we  made  significant  advances  in  adopting  and  utili/.ing  remote  data  for  the 
purposes  of  geologic  and  soils  analysis.  The  primary  technique  employed  was  that  of  image 
ratioing,  a  straightforward  processing  method  which  is  not  only  economical  but  which  has  an 
image  format  similar  to  that  of  an  aerial  photo;  in  the  future,  the  process  can  even  be  aceom 
plished  in  real  time. 


By  intioing  pairs  of  mult ispectral  images,  the  following  geologic  features  can  be  discrimi 
mated: 


0)  exposed  rock  areas 

(2)  separate  rock  groups,  including  Pr<  Cambrian  quartzite;  Precambrian  phyllite.  slate, 
and  schist;  basal  quartzite  and  sandstone;  and  intrusives  such  as  rhyolite,  quartz  latite, 
quartz,  monsonite,  and  phonolite 

(3)  large-scale  structures;  faults,  folds,  intrusive  bodies,  and  regional  trends. 

In  addition,  soil  color,  terrain  patterns,  natural  drainage,  and  soil  texture  are  enhanced  by 
ratio  images. 

To  help  make  the  data  processing  automatic,  a  gating  logic  has  been  developed  for  ratioed 
images.  The  basis  of  this  logic  is  to  use  a  known  ratio  range  for  each  rock  type  in  the  prepara¬ 
tion  of  recognition  maps.  All  ratio  values  which  fall  into  a  preset  range  for  a  given  rock  type 
are  selected  and  automatically  placed  in  the  map  location  of  the  rock's  appearance  in  the  scanned 
scene.  Statistical  procedures  based  on  laboratory  measurements  can  provide  a  firm  basis  for 
discrimination  of  rock  types. 
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Those  advances  in  remote  sensing  analysis  are  supplemental  to  the  thermal  and  mineralogical 
determinations  accomplished  earlier  and  described  in  the  previous  annual  report. 

The  techniques  described  in  detail  in  this  report  now  need  to  be  extender  as  follows:  (1)  by 
application  to  more  rock,  soil,  and  vegetation  types,  and  (2)  through  the  simple  addition  of  air¬ 
borne  ratio  logic  in  order  to  obtain  real-time  ratioed  data.  Because  the  procedures  we  used 
must  be  considered  investigative,  each  step  requiring  the  judgment  of  an  experienced  researcher, 
our  first  recommendation  is  that  enough  analyses  be  performed  to  establish  a  solid  experimental 
and  theoretical  basis  for  the  operational  use  of  ratioing.  Our  further  recommendation  is  that  a 
straightforward  airborne  operational  system  be  developed  and  flown  to  serve  potential  users. 
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